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带 跳 的 随机 微分 方程 近似 解 的 收敛 速率 


摘要 


本 文 给 出 了 三 种 求解 带 跳 的 随机 微分 方程 的 数值 解法 。 首 先 针对 一 类 带 跳 的 型 
随机 微分 方程 ， 给 出 了 基于 Euler-Maruyama 法 的 Split-step 算法 ， 在 方程 的 系数 满足 
Lipschitz 条 件 和 线性 增长 条 件 ， 且 时 兆 函 数 满足 某 种 连续 的 条 件 下 ， 证 明了 Split-step 


算法 的 收敛 速率 为 wAy 人 其 次 , 对 不 含有 时 洁 项 的 随机 微分 方程 的 建立 了 SORE, 
证 明了 当 方 程 的 系数 满足 Lipschitz 条 件 和 线性 增长 条 件 时 ，SS9 算 法 近似 解 的 收敛 速率 
BS. Ba, KSO RIA Bl TN AREAL AZ. UEBA T er HRY eG REAL FE 


SAAD FEE IR SE EE. 最 后 , 应 用 kô- Taylor 展开 公式 对 带 跳 随机 微分 


方程 作 二 阶 展开 ,构建 了 Split-step 一 阶 近似 算法 ， 其 收敛 速率 为 1。 对 上 述 算法 都 给 出 
了 一 些 数 值 实例 来 验证 算法 的 有 效 性 。 


关键 词 ，Split-step 算法 ，SS9 算 法 ，116 一 Taylor 展开 式 ， 全 局 误差 ， 泊 松 跳 
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Convergence rate of numerical solutions to SDEs with jumps 


Abstract 


Three numerical methods for solving a stochastic differential equation with jump are 
given in this paper. Based on the Euler-Maruyama method, we introduce the Split-step 
method to a class of the stochastic differential equation with jump. We prove that it is 
convergence under the Lipschitz condition, the linear growth condition and some continuous 


conditions of the delay function. The convergence rate of Split-step numerical solutions to the 
stochastic differential equation with jump has order & ayaz: Then we discuss the SS8 
method of the stochastic differential equation without delay. We also show that its 
convergence rate has order 5 under the Lipschitz condition and the linear growth condition. 
Next, we expand the conclusion to the stochastic differential equation with delay. We obtain 


that it is convergence rate 了 of SSO numerical solutions to the stochastic delay differential 


equation with jump. Finally, using /to—Taylor expansion, we obtain the first-order 


approximate numerical solutions of the stochastic differential equation with jump. And the 
rate of convergence is 1. Numerical experiments are simulated to testify the performance and 
the effectiveness of the above method. 


Keywords: Split-step method, SS@ method, Itô— Taylor expansion, global error, Poisson 


jump 
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第 1 章 绪论 


三 百 多 年 前 , 牛顿 (Newton，1642-1727) 和 莱 布 尼 兹 (Leibniz，1646-1716) 创 立 了 确定 
性 的 微 积分 ， 对 机 械 动力 学 、 天 文学 等 领域 的 物理 现象 进行 定量 分 析 和 研究 。 随 后 ， 微 
分 方程 被 广泛 用 于 刻画 事物 发 展 的 运动 规律 。 然 而 在 诸多 实际 物理 问题 中 ,依据 一 个 物 
理 问题 来 构造 数学 模型 ， 即 把 一 个 物理 问题 转换 成 一 组 数学 方程 ,但 是 有 时 可 能 无 法 用 
一 个 确定 性 方程 来 描述 物理 问题 。 这 是 因为 物理 问题 的 随机 性 和 复杂 性 不 可 避免 地 给 数 
学 模型 带 来 一 些 不 确定 性 因素 。 这 种 影响 可 以 是 来 自 系 统 本 身 的 不 确定 性 ， 也 可 以 来 自 
系统 所 处 环境 的 不 确定 性 。 例 如 在 信和 号 传输 过 程 中 ， 信 道 存 在 噪声 ; 现代 资本 市 场 理论 
认为 股票 期 货 价 格 具 有 随机 性 等 等 ， 所 以 随机 因素 的 影响 是 存在 的 。1902 年 ，Gibbs 在 
研究 统计 力学 问题 过 程 中 首次 提出 了 随机 微分 方程 模型 。 

1918 年 ，Wiener 对 Brown 运动 在 理论 上 作出 了 精确 的 数学 描述 ， 从 而 建立 了 用 随 
机 过 程 语言 描述 的 严格 的 数学 模型 , 通过 对 Brown 运动 轨道 的 性 质 深入 研究 , 提出 了 在 
布朗 运动 空间 上 的 测度 和 积分 , 这 些 工作 使 得 布朗 运动 及 其 研究 得 到 了 迅速 而 深入 的 发 
展 。 由 于 Brown 运动 也 用 来 专 指 Wiener 过 程 的 数学 模型 , 因此 Brown 运动 也 称 为 Wiener 
过 程 。1923 年 ，Brown 就 Brown 运动 给 出 了 一 系列 严格 的 研究 成 果 。 

1934 年 至 1938 年 S.Bernstein 正式 引入 了 随机 微分 方程 ， 并 证 明了 该 方程 所 确定 的 
随机 变量 的 极限 分 布 是 Kolmogorov 方程 的 基本 解 。 随 后 Gihman Æ S. Bernstein 的 框架 
下 ， 建 立 了 随机 微分 方程 的 相关 理论 。1942 年 It6 在 常 微分 方程 的 基础 上 加 上 一 个 随机 
于 扰 项 ( 即 Wiener 过 程 ) 讨 论 了 如 下 形式 的 随机 微分 方程 

dX (t)= f(t, X (t))dt+ g(t, X (t))dW,, 1>0 
~ =X), 
通过 对 随机 微分 方程 的 研究 得 到 了 关于 Markov 过 程 的 Kolmogorov 方程 ，1944 年 建立 
了 与 常 微分 理论 中 的 Riemann—Stieltjes 积分 不 同 的 It6 积分 ， 即 关于 Wiener 过 程 的 随 
机 积分 , 并 于 1951 年 ,发 表 了 具有 划时代 意义 的 著作 《On stochastic differential equations》。 
随机 积分 为 随机 微分 方程 的 发 展开 辟 了 新 的 天 地 , 而 且 为 随机 分 析 这 门 数 学 新 分 支 的 创 
立 和 发 展 葛 定 了 基础 ， 这 个 理论 被 誉 为 “随机 王国 中 的 牛顿 定律 ”。 

之 所 以 说 It6 积分 为 随机 分 析 的 创立 和 发 展 的 基础 ， 是 因为 对 于 更 接近 物理 直观 的 
郎 之 万 方程 
dX, =E gy +adt , 

m 
这 类 的 扩散 方程 是 很 有 吸引 力 的 。 可 惜 的 是 布朗 运动 WW 对 几乎 所 有 的 轨道 处 处 不 可 微 ， 
所 以 从 数学 上 来 讲 aW 不 能 按照 一 般 的 想法 定义 积分 。 但 是 自从 伊藤 清 给 出 了 积分 的 定 


义 和 随 机 积分 (微分 ) 方 程 后 ， 可 以 用 随机 积分 (微分 ) 方 程 
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X, =X, + [als,X,)ds+ [ b(s,X,)aW, 


的 解 来 表示 扩散 过 程 。 至 今 ， 由 伊 芯 清 开创 的 随机 微 积 分 不 仅 适 用 于 扩散 过 程 ， 而 且 适 
用 于 一 大 类 非常 广泛 的 随机 过 程 半 著 ， 成 为 随机 分 析 的 有 力 工具 。 

可 以 说 ， 随 机 微分 方程 在 金融 系统 、 数 量 经 济 、 控 制 系统 、 统 计 物 理 、 系 统 生物 学 
中 都 有 着 非常 重要 的 作用 中 人。 例如 ， 用 随机 微分 方程 来 解决 期 权 定价 问题 。1973 F, 
Fischer Black 和 Myron Scholes 利用 无 风险 投资 理论 和 随机 微分 方程 理论 ， 得 到 了 著名 
的 Black-Scholes 随机 微分 方程 : 

dX (t)= udt + odW,, 

1 (0) =X, 


并 利用 相应 的 边界 条 件 和 概率 方法 得 到 了 欧式 看 涨 ( 跌 ) 期 权 价格 的 计算 公式 ， 开 拓 了 
金融 工程 从 定性 分 析 进 入 定量 分 析 的 时 代 。 但 是 由 于 随机 系统 本 身 的 复杂 性 ,一 般 情况 
下 很 难得 到 方程 理论 解 的 显示 表达 式 , 只 有 一 些 特殊 的 随机 微分 方程 才能 求 出 其 解析 解 ， 
因此 人 们 只 能 通过 数值 方法 来 研究 方程 的 解 及 其 性 质 。 故 对 随机 微分 方程 数值 求解 方法 
的 研究 是 有 意义 的 。 衡 量 一 个 数值 解法 的 好 坏 主要 基于 三 个 方面 的 考虑 : 〈 一 ) 数值 解 
法 的 收敛 速度 ， 即 数值 解 逼 近 真 解 的 快慢 ; (二 ) 数值 解法 的 计算 复杂 度 ;〈 三 ) 数值 解 
法 是 不 是 能 够 很 好 保持 原 系统 固有 的 性 质 ， 诸 如 稳定 性 ， 散 锡 性 ( 指 系统 具有 一 个 有 界 
吸引 集 , 从 任意 初始 条 件 出 发 的 解 经 过 有 限时 间 后 进入 该 吸引 集 并 随后 保持 在 里 面 )， 守 
恒 性 质 ， 随 机 情形 的 不 变 测度 等 等 。 

随 着 数学 家 们 的 不 懈 努 力 和 计算 机 计算 能 力 的 快 随 提 升 , 对 随机 微分 方程 的 数值 求 
ESET RAMEE TO), REF It6 型 随机 微分 方程 的 数值 求解 有 许多 研究 ， 按 算法 
分 类 主要 有 : 

© Euler 方法 Burrageli1 等 人 讨论 了 处 理 刚性 随机 系统 的 复合 Euler 方法 ， 证 明了 
该 方法 对 于 非 线性 方程 同样 有 效 。ZhangP9 等 人 讨论 了 线性 随机 微分 延迟 方程 的 裂 步 倒 
向 Euler 法 ， 证 明了 该 方法 的 收敛 阶 数 为 1/2 ， 给 出 了 不 同情 况 下 对 于 步 长 取 值 的 限制 。 
Gardón mI 讨论 了 对 于 给 定 的 y={0.5，1，1.5 …} 带 跳 的 随机 微分 方程 Euler 法 的 
随机 Taylor 展开 式 的 1 阶 显 式 ， 并 证 明了 强 Taylor 近似 解 的 收敛 阶 数 为 1。FengP29 等 人 
讨论 了 带 跳 随机 微分 延迟 方程 的 强 Taylor 展开 法 。WangP 等 人 讨论 了 带 跳 随机 微分 延 
IRF EWE Bask Euler 法 ， 证 明了 该 算法 的 收敛 性 和 稳定 性 。 

@ Euler-Maruyama 法 Wangl1 等 人 讨论 了 正 向 的 裂 步 Euler-Maruyama 法 和 裂 步 


Milstein 法 , 并 证 明了 改进 后 的 方法 比 原始 方法 更 稳定 ,Liu 和 Li?" Kubilius 和 Platen 
PUY Protter 的 等 人 研究 了 Euler-Maruyama 法 的 弱 收 敛 性 。Bao 站 等 人 证 明了 对 于 


Vp>2，Euler-Maruyama 法 数值 解 的 p 阶 收敛 速率 为 1p 。Bao 加 等 人 讨论 了 一 类 随机 
延迟 微分 方程 的 Euler-Maruyama 法 ,此 类 方程 的 时 清 项 系数 可 以 不 满足 线性 增长 条 件 ， 
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且 证 明了 该 算法 在 布朗 运动 情形 下 的 收敛 阶 数 为 /2, 纯 跳 情形 下 的 均 方 收敛 阶 数 为 1/2。 
@ Milstein 法 Omar03 等 人 讨论 了 随机 微分 方程 的 复合 Milstein 法 ， 通 过 一 个 指标 
函数 来 决定 计算 中 的 每 一 步 是 采用 半 隐 式 Milstein 法 或 隐 式 Milstein 法 。Higham 中 通过 
作 图 的 方法 给 出 了 线性 方程 的 随机 8 方法 和 半 隐 式 Milstein EMT eK. SKINS 
(3) 考虑 随机 微分 方程 Milstein 方法 的 几乎 必然 及 矩 指数 稳定 性 ,给 出 了 当 步 长 趋 于 零 时 
极限 意义 下 随机 微分 方程 Milstein 方法 的 稳定 性 ,并 证 明了 在 一 定 条 件 下 显 式 和 半 隐 式 
Milstein 方法 都 具有 这 些 稳 定性 。 李 启 勇 和 研究 了 随机 微分 方程 改进 的 半 隐 Milstein 方 
法 的 均 方 稳定 和 渐 近 稳定 性 。 对 线性 检验 方程 ， 得 到 了 改进 的 半 隐 Milstein 方法 均 方 稳 


定 的 充 要 条 件 是 <9< 1。 对 任意 步 长 A> 0， 证 明了 当 方法 的 步 长 充分 小 时 ， 方 法 能 


保持 原 系统 的 渐 近 稳定 性 。 郭 谦 59 给 出 一 个 新 的 求解 线性 随机 时 灌 微 分 方程 的 显 式 分 
裂 步 长 Milstein 格式 ， 运 用 Ito-Taylor 展开 式 证 明 该 格式 相对 于 已 有 的 求解 随机 时 沸 微 
分 方程 的 分 裂 步 长 方法 而 言 具有 更 好 的 收敛 性 。 

@ 0 方法 Ding04 等 人 讨论 了 基于 Euler -Maruyama 法 的 随机 4 方法 , 并 证 明了 该 方 
法 的 收敛 阶 数 为 /2 以 及 稳定 条 件 。Safiquet "等 人 讨论 了 刚性 随机 微分 方程 的 全 隐 式 随 


机 68 法， 通过 参数 9 来 调节 该 方法 的 稳定 区 域 ， 并 且 通 过 数值 算 例 阐释 了 不 管 是 对 于 普 
通 的 线性 方程 还 是 化 学 郎 之 万 方程 亦 或 Duffing-van der Pol oscillator 模型 ,全 隐 式 随机 0 
法 都 优 于 现存 的 方法 。Wange 等 人 讨论 了 带 跳 的 随机 微分 方程 的 补偿 9 方法， 通过 变 
换 将 6 方法 扩展 到 了 带 跳 的 情况 ， 证 明了 该 方法 的 收敛 阶 数 为 1/2， 同时 也 考虑 了 算法 


的 均 方 稳定 的 性 质 。 周 立 群 9 等 人 研究 了 线性 随机 延迟 微分 方程 复合 9 方法 的 收敛 性 。 

© 随机 Runge-Kutta 方 法 Wang! hY T Stratonovich 随机 微分 方程 的 三 阶 强 1 阶 
的 随机 Runge-Kutta 方法 ， 并 具体 写 出 了 一 个 显 式 三 级 随机 Runge-Kutta 方法 和 一 个 半 
隐 式 三 级 随机 Runge-Kutta 方法 。 针 对 It6 型 随机 微分 方程 , Tocino 和 Vigo-Aguiarf 给 
出 了 随机 Runge-Kutta 方法 弱 2 阶 条 件 。2006 年 ，R5gBlerEa 给 出 了 针对 It6 型 随机 微分 
方程 随机 Runge-Kutta 方法 弱 意义 下 的 系数 条 件 ， 并 且 显 式 得 到 了 一 些 新 的 2 阶 随 机 
Runge-Kutta 方法 。2009 Æ, Debrabant 和 R6BlerB4 给 出 了 一 类 更 为 有 效 的 2 阶 随 机 
Runge-Kutta 方法 。 

© Taylor 方法 1982 年 Platen 和 Wagner Xt REN. Taylor 方法 的 各 种 收敛 性 给 出 了 
理论 研究 ，1992 年 Kleoden 和 Platen! KHE It6 公式 给 出 了 随机 Taylor 展开 式 的 构 
造 过 程 ,同时 指出 了 随机 Euler 方 法 可 以 看 成 是 对 随机 Taylor 展开 式 在 0.5 阶 处 的 截断 ， 
Milstein 法 可 以 看 成 是 对 随机 Taylor 展开 式 在 1 阶 处 的 截断 , 而 求解 Stratonovich 型 随机 
微分 方程 的 随机 Runge 一 Kutta 方 法 59 也 是 基于 随机 Talyor 展开 式 而 得 到 的 .2001 年 Tian 
和 BurrageB7 研 究 了 与 Euler 方法 、Milstein 方法 和 1.5 阶 随机 Taylor 方法 相对 应 的 三 种 
隐 式 随机 Taylor 方法 。 


第 4 页 华东 理工 大 学 硕士 学 位 论文 
© Split-step 算法 页 俊 梅 中 ] 给 出 随机 微分 方程 的 Split-step 欧 拉 格 式 的 算法 , 并 证 明 
了 当 方程 的 偏 移 系数 和 扩散 系数 均 满 足 线性 增长 条 件 和 李 普 希 效 条 件 的 情况 下 ， 此 方法 


用 以 求解 随机 微分 方程 的 收敛 性 ， 并 且 求 出 强 收 敛 的 阶 是 /2， 同 时 证 明了 Split-step 近 


似 解 的 均 方 收敛 理论 。Wangt 等 人 给 出 了 自治 标量 Split-step 欧 拉 算法 的 数值 格式 ， 并 
证 明 它 的 收敛 性 和 稳定 性 。 李 启 勇 拓 等 人 研究 一 类 改进 分 步 向 后 Euler 方法 求解 随机 延 


迟 积 分 微分 方程 的 均 方 指 数 稳定 性 。 证 明了 在 约束 网 格 下 ， 该 方法 依 步 长 = 二 保持 原 


系统 的 均 方 指数 稳 定性 。 郝 四 1 讨论 了 求解 随机 延迟 微分 方程 的 分 步 向 前 Euler 方法 在 均 
方 意义 下 的 收敛 性 和 稳定 性 。 将 分 步 向 前 Euler 方法 应 用 于 具有 一 般 形式 的 随机 延迟 微 


分 方程 ， 得 到 差分 格式 ， 证 明 该 格式 在 均 方 意义 下 的 收敛 阶 为 2 ， 给 出 保证 差分 格式 


均 方 稳定 的 步 长 限制 条 件 . 死 园 搬 Gu 讨论 了 带 跳 的 随机 延迟 微分 方程 的 Split-step 算法 ， 
证 明了 该 算法 的 收敛 性 ， 同 时 也 给 出 了 算法 均 方 稳定 的 条 件 。 

@ Split-step 6 算法 (简称 SS 算法 ) Cao!) 等 人 研究 了 随机 延迟 微分 方程 Split-step 
9 算法 的 均 方 稳定 性 和 收敛 性 ， 同 时 证 明了 在 一 定 条 件 下 Split-step 0 算法 的 指数 均 方 稳 
定 并 以 阶 数 0.5 收 你 .Rathinasamyg9 研究 了 线性 随机 延迟 积分 微分 方程 的 Split-step9 算 
法 的 稳定 性 ， 其 中 维 纳 过 程 是 一 个 两 点 分 布 的 离散 型 随机 变量 。 唐 占 涛 Ml 等 人 研究 了 随 
机 分 步 9 CSST) 方法 应 用 于 随机 延迟 微分 方程 (SDDEs) 时 的 稳定 性 性 质 ， 给 出 在 线 
性 增长 条 件 及 单 边 Lipschitz 条 件 下 , SST 数值 解 能 保持 原 方 程 真实 解 儿 乎 必然 指数 稳定 
的 一 个 充分 条 件 。 


ASOT RTE REL (OF (Agr?) Es DHE 116 型 随机 微分 方程 


aX (t)= f(t, X(t))dt + g(t, X O)aW, + h(t, X(t ))aN,, te[0,T] (1-1) 
X(t) = Xo 
或 者 
aX (t)= f (t, X (t), X (t-1))dt+g¢(t,X (t), X (t-1)) aw, 
+h(,X(r),X((t-2) ))an, tefo] (1-2) 
X (t)=(t), te[—r,0] 


LPNS 是 递增 的 ， 右 连续 且 . 和 包含 所 有 的 零 测 集 。X(1) 是 定义 在 [-a.7] 上 的 右 连 左 
极 的 函数 族 ， 其 范 数 | 加 = sup_,。。 


oj， 奢 是 .多 适应 的 布朗 运动 ，N, 是 定义 在 概率 空 
间 上 的 参数 为 4 的 一 维 泊 松 过 程 ， 酌 ，N, 与 .元 相互 独立 。 


AME EER SRE BSH Split-step 算法 的 基础 上 ， 讨 论 方程 (1-2) 的 Split-step 算 
法 近似 解 的 收敛 速率 。 与 裘 半 辰 B11 所 作 的 局 部 误差 不 同 ， 本 文 将 讨论 方程 (1-2) 近 似 解 的 
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全 局 误差 。 在 文献 的 所 研究 的 方程 基础 上 引入 了 泊 松 跳 过 程 , 并 构造 了 Split-step9 方法 ， 
本 文 将 研究 Split-stepO 算法 的 收敛 速率 , 随后 将 结果 推广 到 方程 (1-2)。 最 后 依据 文献 “3 
给 出 得 随机 Taylor 展开 式 和 文献 中 给 出 的 离散 的 随机 Taylor 展开 式 ， 对 方程 (1-1) 构 造 
强 Taylor Split-step 法 。 

本 文 将 如 下 安排 第 二 章 讨论 一 类 带 跳 的 tò 型 随机 微分 方程 的 Split-step 算法 及 其 
收敛 速率 。 第 三 章 给 出 带 跳 的 随机 微分 方程 的 Split-step0 算法 ， 并 证 明 算 法 是 收敛 的 。 
第 四 章 研究 带 时 滞 的 随机 微分 方程 的 Split-step9 算法 。 BABA Taylor 展开 式 建立 一 阶 
的 Split-step 算法 。 
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第 2 章 带 跳 的 变 时 滞 随 机 微分 方程 裂 步 法 近似 解 的 收敛 速率 


本 章节 考虑 如 下 形式 的 具有 时 变 时 滞 的 带 跳 随机 延迟 微分 方程 
dX(t) = f (t, XO), X (t-1(t))) dt + g(t, X ®©), X (t-7@)) aM, 
+h(t, XE) X (0-10 ))av， t e[0,T] (2-1) 
X(t)=9(t), tel[-r,0] 
HF oee Cd-r, 0R) > supa Elo! <o BAFRA, pe) SW, NAA 
Ahi, X(t) ARF C((-r,7];2) MBAS: WAS 适应 的 标准 维 纳 过 程 ， N, 为 参数 
为 大 的 一 维 泊 松 过 程 ， 于 与 W 相互 独立 。 开 (六 ) = lim X (s) 。 方 程 (2-1) 的 时 滞 项 7z(1) 为 
时 变 时 洁 ， 假 设 时 洁 r(f) 满足 下 列 两 个 条 件 ; 
(1) rt): [0o) 一 人 ， 存 在 某 个 正 实 数 >， 使 得 t+-r(D)>-r。 
(2) 对 任意 的 上 s20， 存 在 一 非 负 实 数 p ， 使 得 |r)-zrGsjls plt—s 
为 了 保证 方程 (2-1) 存 在 唯一 解 及 后 续 证 明 需 要 ， 我 们 引入 如 下 假设 : 
假设 2.1 初 值 函 数 g 关 于 指数 y 连续, 即 对 于 任意 的 s, fe[-r,0] ,存在 正常 数 天 ， 
使 得 
El-p) < Kle- s”. (2-2) 


假设 2.2 方程 系数 /，g， 关于 变 元 /满足 指数 4 连续， 关于 变 元 x，y 满足 
Lipschitz 条 件 ， 即 对 任意 的 t,t, Z0URn, n y eR, BEES a>0. Ef 
BK, 44 
|f ys) — fb, x Ya) V g(t, x, 7)— Bla» xX yy hh, ,7) -Alx 

< K,(t,-t,|" +|x, -—x,|+|», -y 

假设 2.3 方程 系数 f ，g , h 均 满足 线性 增长 条 件 , 即 对 任意 的 :>0 以 及 x，, ye RR， 
存在 非 负 常数 天 ， 使 得 
If») vlet viaa K a+? +p. 

由 文献 中 可 知 ， 当 假设 2.2 和 假设 2.3 成 立时 ， 方 程 (2-1) 存 在 唯一 解 。 
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现在 给 出 方程 (2-1) 的 Split-step 算法 。 首 先 ， 我 们 定义 一 簇 在 [0,7] 上 的 标准 网 格 点 
I= <h < < < =T}; HHL =kA, ke(0,1,2,...,N)> aE. 由 于 与 
时 沾 相 关 的 点 t. -z(t) 不 一 定 会 落 在 标准 网 格 点 上 ， 所 以 我 们 针对 所 有 的 s, =1,-7(t,) 
定义 一 簇 非 标准 网 格 点 gpy 会 {-r < sy <s <0 <5, <…<sy ST} © 
将 上 述 两 个 网 格 点 合并 ,得 Ns 全 {~r <5) <5, <…<h <…<sy <…<ty =T} 
显然 ， 对 任意 的 seN;， 当 s, >0 时 ， 必 存在 te S、 和 ce (0, 切 ， 使 得 s, =1,+ ARIZ. 
用 Y(s,) 表示 X(s,) 的 估计 值 ， 对 于 任意 的 seN;， 当 5s,<0 时 ，F(s,)=g(s,); 而 当 
s, >0 时 ， 用 于 计算 方程 (2-1) 数 值 解 的 Split-step 算法 可 表示 为 
Y, =Y(t,) + Aft, Y, Y, -t(,))). (2-3) 
Yll, +A)= Y; + glt Yy Y, —T(t NMG +A) -W t) 
+h(t,, Y, Y(t, 16, DENCE, + A) -NG )). (2-4) 
特别 地 ， 当 s = + 的 落 在 标准 网 格 点 3S 中， 即 @=1，s, = 时 ，Y(&w) 的 计算 


式 为 
Y; = Y(t.) + Af (te Yrs Yt -7(t,))), (2-5) 


¥ (ty) = X: +g(t,; Y > Ter, —T(t, JAW, + h(t,, T ， Y(t, —t(t,)))AN,. (2-6) 


2.1 全 局 误差 的 收敛 速率 


本 节 将 证 明 Split-step 算法 的 收敛 性 ， 即 算法 (2-3) 和 (2-4) 的 数值 解 一 致 收敛 到 方程 
(2-1) 的 精确 解 。 下 面 先 给 出 Split-step 算法 的 全 局 误差 和 均 方 意义 下 收敛 的 定义 。 


定义 2.1 对 所 有 的 s,e Ns， 必 存在 i e Sy 和 ce (0,1]， 使 得 s, = 大 + 弘 ， 算 法 (2-3) 
和 (2-4) 的 全 局 误差 e(s,) 定 义 为 e(s,)=XY(t, +cA)- Y(t, + 6A)。 
定义 2.2 设 久 (1) 为 方程 (2-1) 的 精确 解 ，7 (1) 为 通过 数值 方法 得 到 的 近似 估计 值 , 若 


FETE BRC Hy, BE sup [x ()—¥ (0 |Z |< C-A” RIL, 则 称 在 均 方 意义 下 近似 
tel0,7 


解 Y(t) 的 收敛 ， 其 收敛 速率 为 7 。 
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现在 给 出 本 章 的 主要 结论 。 


定理 2.1 者 方程 (2.1) 满 足 假设 2.1-2.3; MA<minli i i 则 存在 与 步 
1 


J4K, 
长 A 无 关 的 常数 C 时 , 使 得 不 等 式 ， le( sup |aOF x] < CAT"? RAT, 即 通过 算法 (2-3) 


和 (2-4) 得 到 的 数值 解 以 速率 a Ay 人 一 致 收敛 到 方程 (2-1) 的 精确 解 。 
为 了 证 明 上 述 收敛 定理 ， 我 们 需要 引入 下 面 的 几 个 引 理 。 
引 理 2.1 如 果 假 设 2.3 成 立 ， 则 对 任意 的 te[-r,7]， 其 中 0<7T<oo， 存 在 正 数 C ， 


使 得 E( sup XOF A) <Cd+zlof), 
其 中 多 (1) 为 方程 (2-1) 的 解 。 
证 明 ; 由 方程 (2-1) 可 知 ， 对 于 所 有 的 0<+:<T， 有 
XO)= XO+ | fu, X), Xu ra) aut | gu, Xu), Kur Nam, 
最 [ hlu, X (u), X (u -r(u)))dN,. 
两 边 先 在 [0%]j 上 取 上 确 界 再 取 期 望 ， 由 基本 不 等 式 


(a+b+c+d)y <4a°+4b°+4c° +42 可知 ， 对 于 所 有 的 se(0,T]， 有 


E( sp xof F< 4{E|X (Of + E(sup | [ f(u, X (u), Xu—ru))) du |Z) 
Ost x Ost<s, 
2 
+ E(sup|[ (u, X (1), X(u-t(u)))aM,} |F) 
2 
+ Elsup|[ hu, X0), X (u -1AN | FD) 2-7) 


由 柯 西施 瓦 兹 不 等 式 ，Fubini 定理 和 假设 2.3 可 得 ， 


2 


Elsup| (fe X0) Xu- |.) 


< sB F IXe, Xr OA | 


<K,s, jà E[( sup HXP +|X@-c(@))|)| Hau 


i 


< K,5? + Ks, f =| (sup (XOF +|x(@-c@)) 
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< K,s? +2K,s, | E( sup Xo [gan (2-8) 


由 Doob APEA, FALMER ATA AT en, 


2 


E(sup || glu X(w),X(u-ru)))aW, 
Ost<s, 


A) 
<4 Edet, Xu), Xu- | Hau 


<4K,[s, +2 f E( sup |X(@P Ls 
~rsOsu 
= 4K,5, +8K, [' EC sup xo |x (2-9) 
—rsOsu 
对 于 泊 松 跳 , WN ,=N,-d4u, HB-D-GASRB 


2 


A) 


E(sup | (Hu, XW), Xu-e@))dN, 
SS 


2 
= E( = | [hu,¥ (u), X (u —t(u)))dN, +A (aq, X(u), X(u- r(u))) do a 


2 


f hu, X0), X(u-c(u)) aN, 


A) 


2 


<2E(sup 
Ostss, 


+24? (sup | | au, X (U), X (ur) F) 
Ost<s, 


z| + zasa f lau, X (u), X (u — r(u)))| cu) 


<8AE [ r |a(u, X (u),X (u ~r(u)))| du) 


<22K,(4+ Asis, +2 | E( sup |X (oj 


将 式 (2-8)- (2-10) 代 入 式 (2-7)， 有 


E(sup |X(OF 
Ostss, 


Z) < 4{E( sup |p(O)| ) + K,s,0, +2K,0, f ' E( sup XF 
-rsĝ<0 -rs@su 


其 中 O =T+4+2A(4+AT). 
注意 到 

Sup |x wy < ap LOF v mp [x wf ; 
故 由 式 (2-11) 可 得 


F jdu]. (2-10) 


F)du} (2-11) 
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E( sup xof F )du}. 


F) S HE( sup PO) + K,s,0, +2K,0, [ E( sup XOF 
-rsĝs0 =rsĝsu 


由 Gronwall 不 等 式 ， 可 得 


E( sup |X (OF 
—rstsT 


T) SAE, sup |p) )+ K,s,0, 10 + 8") 
-rs800 


<4(K,TO, v D0 +e A1 + Elo). (2-12) 
引 理 2.2 如 果 假 设 2.3 成 立 , 则 对 任意 的 0<s <t<7, 其 中 0<T<w, Ht-s<1, 
存在 正 数 C, ， 使 得 方程 (2-1) 的 解 满足 妃 (|x(D-XCGsjP| 画 jsCC-9。 


证 明 : 由 方程 (2-1) 和 基本 不 等 式 (a+b+c) <3a +3b°4+3c?, WH 


s) 


E(|XO-X¢f 


4) «ss| [7G xe, Xu -rad 


s) <3K,(t-s)l szmal] i +X +X (u -r(u))} Ju 


s| 
下 | 


2 


:al fe, XC), Xu rl)am, 


2 


+ | [ nu, X00), Xu t(DaN, 


用 与 引 理 2.1 类 似 的 证 明 方 法 可 得 


E(IKO-XGjf 


+3K, (+2aye| [itlx(f +|X(u-2(w)f au) 
由 引 理 2.1 UK t-s <1 可 知 


E(|X@-X(f 


x) <6K,(1+4+ ai + 26 sup XOF 
—rstsT 
s C, (t i s), 


其 中 C =6K,0+4+4 1+ 2C,1+ Elo). 

为 了 证 明定 理 2.1。 RIK X N) ， 和 (一 z(G)) 分 别 代入 算法 (2-3) 和 (2-4) 的 右边 的 
F, Y-t) RITIKA: 
Y; =X (t) +A (tY; X(t -7(t))), (2-13) 


Yit, +Sh)= Y; + a(t ¥, X (t, -Th WG, +A) -WA 
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+All, Y, X(t - T(t, ING, +A)-NG)I- (2-14) 
对 (2-13) 式 中 的 立 作 估 计 。 
引 理 23 如 果 假 设 23 成 立 ， 则 对 任意 的 ze[07]， 其 中 0<T<o ， 当 


amin az} FELETEM LC,» Co (AESERE(2-13)-P BRE DLN Y, WE 


a %)sCElof +C.. (2-15) 


i 
al 


2 
+ 
r 


证 明 : 对 算法 (2-13) 式 两 边 平方 后 取 期 望 ， 由 基本 不 等 式 (a+b) <2a? +2b? 可 知 


|S) +2867? (P(x (4-7 (%))) 5). 


K a) < 28 (|x¢, 
由 线性 增长 条 件 可 知 


a) 


rp 


F)+E(|xG —t(t, 


valk 


I) <2 (XG |.%)+26°"K, (1+al 


移 项 整理 后 ， 由 引 理 2.1 可 得 
2 2 
5) < ry (x )| |5) 
26°A°K. 2 
Tota (RG reel] %)+ 
<20+67A°K;) 
1-27A'K, 
<C,Elof +C, 


26°A’K, 
1-26°A’K, 
26°A’K, 


2 
(1+ Eljo| )+122A5K 


其 中 C; =4(1+ K,),C, =C,+4K,. 


引 理 2.4 BX (t) 272-1) HH, YO Split-step 算法 得 到 的 近似 解 ， 若 假设 


ee a a S 
2.1-2.3 RUIZ, MAER OSksN-1, 存在 正常 数 C,， A<min| et 


时 ， 
E(le@, + 5A) A) SCA. 
证 明 : 利用 方程 (2-1) 和 式 (2-4)， 得 到 误差 
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E(t, +A) = X(t, +cA)— Y(t, +A) 
=e(t,)+ [CXX E-D- 16, Ty 2G, -rE lar 


+ [Le XO), X-7O)- gl, Y, Yit, -16 IAM, 
+ [Ar XO), Xr rr) h(t, ,Tt -T(t AN, (2-16) 
HRN +A 运用 It6 公式 得 


g(r, X (r), X (r -10))- (4. TI -ro dr 


E (t, +A) = 67 (t,)+ me 
+ “hr, XOX E-t) -Alt Yy Yl, —t(t, Di dN, 
+2 ie Lf. X(r),X(r—-c0r)) Sto T Ilt, -7(4,) Near 
+2 f [g(r, X(r), X(r-1(r) - g(t, T I, - 7, NeW, 
+2 i [h(r, X) X (r-r) — h(t, T; Vt, — 7, le(AN, (2-17) 
再 由 E( [soaw]|%)=0, E(f/@an,|.%)=28( [roa] %), 2-1 PHA 


条 件 期 望 得 
E(t, + A)%) = E(P) 


+B| (C lee XOXE -r-e T; 20, -aD ar] | 


5 


+ | f Ahr, XE), XE -10 -hlt ,TC -10| dr 
*| 


| 


+22 [Lexo -r(r)- f(t, Ys YG, - 7, le) ar 


+2AE | ine (iht, X(r),X(r —t(r)))— hlt, Y, , Y(t -—t(h e(r)) ar 
(2-18) 
现在 估计 (2-18) 式 的 右边 第 二 项 ， 由 Lipschitz 条 件 和 基本 不 等 式 可 得 


allay 3) 


< 22( [OXOX E-t- 86, %5.XG- eG] ar 


g(r, X(r),X(r-2(r) g(t, Ty Tl, -T(t dr 


z) 


s) 


3 z( 4 BORA ,Xt -Tt -g(t, ¥, , Y(t, -ry ri 


3 2 
<2K7E 上 dr 


o 


lr 1," + Xe) -Y; |+|X(r -r(r))- X(t, -r(t,))| 
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5 


2 tA 
+2K7E * 


= 2 
-r(t,))- Yit, -rt dr 


<6K? [ie 人 ra elo- Eg 


+6K? [a(x (r—c(r)) - X(t, -TED IA 


+4K? ls f Ha +E|| -rcD)-Ze -rf | jar (2-19) 
首先 信 计 (| I)e ATER (a+b) <2a? +26". IE 22 和 
式 (2-13) 可 知 


| 


Jaja- 人 zx-xzcD+xeo-z 


“| x) 


sof" hes: x(t, )P Ajar” z| | 四“ 
<2C, i “b= |ar+2 f (py GY XG - z )a. 
利用 假设 2.3、 引 理 2.1 和 引 理 2.3 可 知 
| sfl) 
sp F}+E(|XG,- c(t, |z) 
<C,A? +2K,A°[C, Elof +C, +C,0+ Ello) +1) 
<C A +A . (2-20) 


其 中 G =2K G +C, Ell +C, +C, +]. 


— F 


j 5); 由 算法 (2-3)、(2-4)、(2-13)、(2-14) 和 基本 不 等 式 可 知 


a(x- 


F) =E (XG) + NG XG -D 


5) 


+ 去 了 站 2 
< B( {ot )+K, Aly; -Y, a ok ii 


ig 


AFE) + oAf tE, YG -DY 


s). 


<3E(|e()f 


%)+3K; ez 人 


+3K7°E (le, -zt ‘| 5) 


移 项 整理 得 
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3K; A’ 


如 i TS 


s) szga Plet if |%}+ E(le, -eDi (2-21) 


接 下 来 需要 估计 【『 ”2(|xe -rc)- XG, -a 
its 
(1) #r—z(r)<t, —r(t,) <0 RY -7(,) sr—c(r) <0, WR 2.1 可 知 


F )dr ， 分 以 下 三 种 情况 进行 讨 


zl 人 lz -t(r)-X(t, -T(t DFL) <K(r-70)-@ -7@ ) 
< K({\r-4,|+|-@)-7@,)[P") 
< Kir —t,|+ plr -tb 
即 
E( [Xr -rr) - xX, rt) L5) <K(+p) A”. (2-22) 
DE O<r—r(r) <t, —1(t,) MO St, 一 rt(f)<r 一 tr(r)， 则 由 引 理 2.2 可 知 
E(|X(r-r0) -Xr I A) < C(e- -t 76) 
即 用 上 面 同样 的 方法 可 得 
E\X(r-2(r))- X(t, -rt $C, 0+ A. (2-23) 
BË r—t(r) <0 <t, —r(t,) Rt, -rlt <0 <r-t(r) BEB Er-z(r)< o< t, -T(t,), 
HR -(r —2(r)) -r-r -1t $ (14+ PAU 
t, -t(t,) <t, -r(t,)-(r-t(r)) $ (14+ p)A RM, BW 
E(|K@-r)-X-7@)|H) 
< 2E(|X¢r -r(r))-X(0)P |5) $ 2 (|X(0) -X(t -r(t 六 | 气 ) 
<2E(|p(r -e-f |A) +(x -XG -eD A) 
<2K|r—r(r)|” +2C, |t, -z(t,)| 
<2K(1+ p)” A” +2C, (1+ p)A. (2-24) 
同 理 ， 当 一 z(t,)<0<r-zt(r) 时 可 以 得 到 相同 的 结论 。 
综合 式 (2-22)、 式 (2-23) 与 式 (2-24), 可 得 
Ef" (K@-2@)-X@ -re LA) drs CAL, (2-25) 
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其 中 C =2(1+ p)” K +2(1+p)C, o 

将 式 (2-20)、 式 (2-21) 与 式 (2-25) 代 入 (2-19) 得 到 
z( [ -T(r)) -g(t, ,Yh -TCh yf dr 


9 


< eE q 46, A? ÓA CA] 
3K; A? 
)+ 2 人 2 
1-3KA 


(le -tN 


nrag rle z) 


s)|. (2-26) 
类 似 地 ， 估 计 (2-18) 式 的 右边 第 三 项 ， 得 到 
| [ he, X(r), X(r-t(r)))-Alt,, Ý p , Y(t, —T(t, Ma 


+E(leG, -r 


< skia An < T+ CA +CA' +C eas] 


+ arial a E (lett, ) 一 二 (ed -r(t 洲 辐 | (2-27) 


1 
1-—3K?A? z) 1-3K2A? 
现在 估计 (2-18) 式 的 右边 第 四 项 ， 根 据 Holder 不 等 式 ， 用 前 面 同样 的 方法 可 得 
| Pi (r, XO) Xr -1r - FG, ¥, Yit, -1i mee 5) 


int 


<7 J E(|F0.X0),X0 -10))-St ,70 


a 


E(le@f 同好 


BIFE XOXE- FV, MG -EDS |) or 
k [zl(eoF ls)s 


Aza! S . ns 
S ala +C,A? +E,A? +C,A% | 
neral 3K?A 2 E(le@,)) |%)+ ry E(leG -7@)f 四 | 
1 
ig fl E(le@ A )ar (2-28) 


类 似 地 对 (2-18) 式 的 右边 第 五 项 作 估 计 ， 可 得 
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225( [ (cx Xr -rmi ,TG -rE le) dr z) 


2a+1l 
< 人 +C,A? +GAs Egor | 
2a+1 
242 3 2 1 2 
+16K; A Al a rle IS) + sm let ~r(t,))| )| 


1 as 
站 二- 


E(le@)) |%) ar (2-29) 
ey 式 (2-27)、 式 (2-28) 和 式 (2-29) 代 入 式 (2-18)， 可 得 
E(|e@, +A) |%) 


< E(lett, f jajar sasl E +CA2+CA3+ | 


+(442 449 44| 开 


和 E (ec = Thi, yf EA )| 


E(le@f |%)+ mer 
1 
二 f E(f L5) dr | (2-30) 
XF s, =t -1(f,)ENy ， 必 存在 整数 <k 且 满足 1: e 3y Ace (0,1), HR 
si =t +h 成立。 于 是 有 E(|e(t -16D A) = sll, +A) |5). 


E(|elt +A 5) 


< Efe |%)+ (44? +4+5)- KE 一 +C A? +CIA3 BY Aerie 


HAAR +4+5)- KA -3pm 


le, rafia) J Elle |x) ar 


ject L5) 


1 
十 一 一 一 
1-3K7A’ cae ( 


< El(leGf |%)+ (44? +44+5)- Kia at C,A? +ÕA +Č AC] 


+ (Ad? +445)-4K Ol pars me Efe tIS) 


"E mee raga (a? (lee, reall 1% pes Ele) |%) ar. (2-31) 


注意 到 
5 fello ias) [°2(XO-X6 + A) + XG, A-F, [Har 


< i E(\X)-XG + A)f |%)ar+ “E(X +A) -F +a) |%) ar 
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< OA +E (let, +cA)P |%)4 (2-32) 


代入 (2-31)， 得 
E(leG + eA) |%) 


2a+1 


S/S ce(0] 2a+1 


2 i 2 3 | 2 ) 
+(44' +4 +5) AKPA — a rages =e lee, +A |% 
TE i Ln : 
+I E: sp ë (let +) 1% Jic +AE (|en +A |-%) (2-33) 


式 (2-33) 移 项 整理 后 ， 左 边 先 去 上 确 界 再 取 最 大 值 ， 令 
E, = max sup E (let ASIA) Fe 
0,1 


Osisk-1 ç 
[1-3K2A7][1-A] 1-A 


zi [ 1 十 16K; A(4A7 + A + ie Å C. ACDA) (2-34) 


1-A [1-3K?A’][1-A] 1-A 
IBAA 
1 16K7A(447+A+5) 
人 
"SGA [1-3KA7][1-A] im 
Cs Aladar S| 1 16K A(4A7 + 1+5) 


PIA T FE J 
1-A 名 1-A [D-3K?A N-A] 
< Cs 和 CarDA2A(274 y 1 Š 16K? A(4A? +4+5) j 
= 气 1-A [1-3K?A?)}[1-A] 
a 16K; A(4A7 +A + S)w A 
< Ce 人 (za+NA2A(C27+1) 1 一 人 fl —3K? A’] fl = A] 
1-A A 16K7A(44°+A+5) 
1-A [1-3K/(A)’][1-4] 
16K7A(4A? +A+4+5 1 
eR Gal 
< C, 人 (zc+lD)^2A(27+1) 1 
16K7A(4/7 +4+5) 


人 十 
1-3K/;A’ 
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16K2A(412 4+S) 
in ee = 
<C 人 ?2aAlA27 [+ y= 3K, PA? ] ir ab 
ie i ea 
pe S 
1-3K2A 
HE A=—, 所 以 当 人 一 0 ， 即 和 一 oo 时， 
2 2 N 2 2 
hessa stes] |) EE S, +] BREE WEF. XM, 


1 caw 
Lal ERR. 所 以 


Elle, +A) 5) < am ; 


引 理 2.5 如 果 假 设 2.1-2.3 成 立 ， 则 对 任意 的 te[0,T]， HHO<T<0o, 4 


1 
sm- 支 | 存在 正常 Bi CC > 使 得 


证 明 : 由 引 理 2.1 和 引 理 2.4 可 知 
z(E [ls =z(leG) -X |S) 
<2E(e(s))f |%)+2E(|X(s)P IE) 
sti. 


| (2-35) 
对 算法 (2-5) 式 两 边 平方 后 取 期 望 ， 由 基本 不 等 式 (ca+ 信 <2a*+2 大 可知 


z 人 区 %)<28 (Pe |%)+2Eea: (f° (6 F, F (t,-7(4.)))|H)- 
由 线性 增长 条 件 可 知 

lal 5)<2E(FG) 5) +20°K, (ual Si 
移 项 整理 后 ， 由 引 理 2.1 可 得 

z 人 zf 列 = T (ref nk A |a) 2 


1-2A°K, 
2(1+ A’K,) 2A°K. 
~ 1-24°K, ” 1-2A?K, 


2 
y 
r, 


5} EPG -re 


z) 2A2K, 
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oC. (2-36) 
其 中 Cs =4(1+ K,)C,+2K,. 
为 证 明 全 局 误差 的 收敛 速率 ， 定 义 连 续 时 间 近 似 解 ， 当 ie [和 ,sa) T EX 
Y() = Y(t,) +(t-4,) St Yi YG, -7G,))). 
+ g(ti, Yp Yt —7(t,))) (AW, - AW, ) + hlt, Yi Y(t, —7(G,))) (AN, -AN,)- (2-37) 


将 (2-3) 和 (2-6) 逐 次 代入 (2-37) 可 得 : 
¥()=%+(t-t,) S tY, Yt, -1t )) 
+8(t,, Fri -T(t (AW, — AW, )+ h(t,» FFG, = t(t,)))(AN, = AN, ). 


+S Nf (0, ,FO tl) + SF 9.7.70, -z(t AW, 
i=0 i=0 
+5 ps ,FG -rl DAN, (2-38) 
接着 定义 3 个 示 性 函数 
NI, _ 
4 (t) = 2, Y, Kin tur) (t) + YyX er) (t), 
Z, (1) 2, Y (t) Zins) (1) Fy Xuan} (1), 
S N- 
i D Rina) (t). 


tg F 


其 中 z= ， 5 。 于 是 式 (2.38) 可 改写 为 


¥(t)=H+ [lez (s).Z, (s —r(s))\ds+ [ez (s),2, (s-z(s)) aw, + Je Z,(s).2; (s- r(s)) AN, 


(2-39) 
接 下 来 ， 给 出 定理 2.1 的 证 明 。 
定理 2.1 的 证 明 : 由 方程 (2-1) 和 式 (2-39) 可 知 


e(t)= | f(s,X(s),X(s—7(s))) -7 (5,2,(5),2,(s -7(s))) es 
$ fel xOx 6-9) -se(52 (s),Z, (s -7t(s))) am. 


+ [a(s,X(s), xX(Gs -7(s) ))-h(5,Z, (s),Z,(s—7(s)) aN, 
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两 边 平方 后 ， 先 取 上 确 界 ， 再 取 条 件 期 望 ， 得 
g 


el sup (af ry < s(a |/ (x06),x(s-7(0)-/ (5 Z,(s),Z,(s -r(s)))}as 
2 1 


(2-40) 
接 下 来 ， 对 于 (2-40) 右 边 的 三 项 分 别 做 估计 。 对 于 te[,)， 根据 Holder FER, 
假设 2.2 可 知 
3 


< sre | fl Ff (s,X(s),X (s—t(s)))-F (5z (s), Z, (s-z (s))) i 


2 


+3E (æ 
Osust 


he (s,X(s),X(s—r(s))) -g(5,2, (s),2,(s—r( 5s))am, 


[a(s,X(s), xX((s—7(s)))-—n(s,2, 


safal 


2 


(Gc X(8), X (s-1(8)))- £ (5,2, (s),Z,(s—2(s))) ds 


3E (æ sup 


9 


zeara fxe- -1(s))-Z,(s-e(s))f ds 


3 


HFS [(s-1)as+ 全 -ae = 和 and? s3NA =37A Ee ae 


t 4, 2 


x 67 K CANA Brean, 


< 37K? (| fis- ses} esx [zeo- Z,(s) ds 


< 31K? s I(e-n)as + Is -t, p)» meS] J \x(s)—Z,(s)f ds + fx(s)-2 (sj ds 
+ are fj (s—t(s))-Z,(s—t (fal 


和 式 (2-21) 可 得 


Je X(s),X(s- 1(s)))- faz 1(5),Z. a(s- -到 


3 五 | sup 
Osust 


3 


242 
<<T’ K7A+67°K? (C,A+C,A’)+ TARA 
2 


1-3K? A? 
7+3K? A? 
1-3K?A? 


3 


SSP KPA+6T’K? (C,A+C,A”) + AT KOA ra 


(2-41) 
先 利用 B-D-G 不 等 式 得 到 
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ý 2 
Ee A | 
Osust 0 


人 g(s,X(s),X (s—r(s)))-8(8,Z,(s),Z(s-2(s))) eM, 
< oil je X(s),X (s-1(s)))—g(8,Z,(s),Z,(s -r(s)))) as | 


再 利用 假设 2.2， 同 第 一 项 的 方法 可 得 


3E E 
Osust 0 


(e(s.X(),X(s-t(9)))-8(5,.Z,(s),Z,(s-z(s))) 


7+3K; A? 


ape xI2TK C ANA (2-42) 
1 


<6TK?A+24TK? (C,A+Č,4? )+ 


. 类似 地 ， 估 计 (2-40) 右 边 的 第 三 项 ， 得 
{wa Ja(s,X(s- ),X ((s- t(s) ))- h(5,Z :(s),Z,(s—z(s) nya | 


ser sul 他 (s, X(s™ ,X((s- T(s) ))- h(5, Z,(s),Z. 2(s- r(s)) Neal | 
vee sl f(a (s, x ),X ((s- r(s)) ))- h(5, Z,(s),Z,(s- -r(s))) as J 


5 


x(244+8A°T)IKPC A" 


<(244+84°T) felh (s, X6), X (Cs—7(s)) )) -a(s,2,(s),Z(s—r O) 


7+3K;A? 


<(124+44°T)TKPA+2TK? (244+ 82°T)(C,A+C,M) +T Ran 


(2-43) 
将 (2-41)- (2-43) 代 入 (2-40) 整 理 可 得 


a( sup ZON x } < ČA? 
OsisT 
由 于 sup ls = sup le (of v sup le (ef ， 运 用 假设 2.1 得 
—rsisT ~rs1s0 OsisT - 
E ( sup (oF EA ) <E (sup le (o EA ) +E ( sup |e of) 
-tst<T Ost<T -rsis0 


< ANN + KAY 
< C A2ertral 
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2.2 ”数值 模拟 


为 了 说 明 本 文 所 给 出 的 均 方 稳定 条 件 以 及 算法 的 有 效 性 , 本 节 将 对 方程 (2-1) 进 行 数 
值 模拟 ,文中 所 给 出 的 近似 解 和 误差 都 是 100 次 计算 得 出 的 均值 。 

例 2.19) 考虑 如 下 形式 的 方程 
| (t)=[ aX (t)+bX (t-1)]dt+cX (t-1)aW (t)+ DX (t)adN (t),t 20 


2-44 
X (t)=t+1,t e [-1,0] aiá 


对 式 (2-44) 的 系数 取 a= -9,2 =7,c=1,D=0.2,A=1. 精确 解 是 根据 [38] 的 方法 得 到 的 ， 即 


取 步 长 1/1024 的 Euler-Maruyama 算法 的 近似 解 作 为 精确 解 ， 图 2.1 给 出 A=0.01 的 
Split-step 算法 的 近似 解 和 Euler-Maruyama 算法 近似 解 ， 从 图 2.1 可 以 看 出 Split-step 算 
法 的 数值 更 近似 于 精确 解 . 图 2.2 显示 的 是 Split-step 算法 的 近似 解 和 Euler-Maruyama 
算法 近似 解 与 精确 解 之 间 的 误差 ， 由 图 2.2 可 知 ，Split-step 算法 的 近似 解 与 精确 解 之 间 
的 误差 更 小 。 此 例 说 明 Split-step 算法 要 优 于 Euler-Maruyama 算法 。 


2.1 精确 解 Split-step 算法 近似 解 、 和 Euler-Maruyama 算法 近似 解 
Fig. 2.1 Exact Solution, Solution of the Euler-Maruyama scheme and of the Split-step scheme 
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‘ —-- RPK 
"E Ni: ATE Ai fet % Mo o3 ae 

afi i i % F 了 i Mi " F “hy ` : F y ki ` Si ! Ns 

2 EN zla i vah i i pe Mt ‘an At 3 1 SA kag l O 
i Why! wA a nA Wi | t “tes 于 
A: ity I Ww! fi i ! m ly 1 Ea ! \ Phy! th Pe 
3 hy 5 1 ul 1 rhb | it ih } y igi 4 my 
Mit nt ui wi 1! i ri eo 
í | f wi if 1 WN 
a | À! \ | 


| 


图 2.2 离散 解 和 精确 解 点 对 点 之 间 取 对 数 后 的 误差 
Fig. 2.2 Point mean square errors of the Euler scheme and of the Split-step scheme in a logarithmic scale 
例 2.2 考虑 如 下 形式 的 方程 

aX (t)=[-2X (t)+ X (¢-1) | dt+ X(t) dW (t)- X (4) dN (t),t 20 
| X (t)=t+1,t e [-1,0] 
图 2.3 的 精确 解 是 取 步 长 1/1024 的 Euler-Maruyama 算法 的 近似 解 作为 精确 解 ， 从 图 2.3 
可 以 看 出 Split-step 算法 的 数值 更 近似 于 精确 解 。 图 2.4 显示 的 是 Split-step 算法 的 近似 
解 和 Euler-Maruyama 算法 近似 解 与 精确 解 之 间 的 误差 ， 由 图 2.4 可 知 ，Split-step 算法 的 
近似 解 与 精确 解 之 间 的 误差 更 小 。 此 例 同样 说 明 Split-step 算法 要 优 于 Euler-Maruyama 
算法 。 


(2-45) 


2.3 A =0.01 Ff Split-step 算法 、Euler-Maruyama 算法 和 精确 解 模拟 结果 


Fig. 2.3 Exact Solution, Solution of the Euler-Maruyama scheme and of the Split-step scheme with 
A=0.01 
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point to point error 


图 2.4 离散 解 和 精确 解 点 对 点 之 间 取 对 数 后 的 误差 
Fig. 2.4 Point mean square errors of the Euler scheme and of the Split-step scheme in a logarithmic scale 
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第 3 章 ” 带 跳 的 随机 微分 方程 39 法 近似 解 的 收敛 速率 


本 章节 考虑 带 跳 的 随机 微分 方程 : 
AX (t)= f (t, X (t))dt+ g(t, X (t) dW, +h(t,X ())daN, relo7] 
X(0)=9(0) 
EOWA SF 适应 的 一 维 布朗 运动 ， N ,为 参数 为 妈 的 一 维 泊 松 过 程 ， 且 于 与 N 相互 独 


(3-1) 


立 . 对 于 所 有 的 :>0 ，y9(0) SW, M N, HAIM, HEO <0; XNARBF 


C([0,T]; R") 的 随机 变量 。 假 设 f,g,h 为 其 所 有 变 元 的 连续 函数 。 
为 了 保证 方程 (3-1) 存 在 唯一 解 ， 还 需要 引入 如 下 假设 : 
假设 3.1 方程 系数 f,g,h 关 于 t 满 足 Holder 连续 ， 关 于 第 二 个 变 元 满足 Lipschitz 


条 件 ， 即 对 于 任意 的 0<t <7,0<t, <T 以 及 x,yeR， 存 在 正常 数 芒 ， 使 得 
| 1 
es) -leveled -gC lv As Me SK (hy -bl + (3-2) 


假设 3.2 系数 f,g,h 均 满足 线性 增长 条 件 ， 即 对 任意 的 0<1<7T 以 及 xeR， 存 在 
正常 数 K,， 使 得 


FEP vl oF va oF Ka+) (3-3) 


由 文献 69 可 知 ， 在 上 述 假设 条 件 下 ， 方 程 (3-1) 存 在 唯一 解 。 
现在 , 将 对 带 跳 的 非 线性 随机 微分 方程 (3-1) 建 立 SS 算法 。 取 和 A 为 步 长 ， 了 = NA， 


其 中 N 均 为 正 整 数 .节点 4 <kA, k=0,1..,N 4-08, 二 J(0); 而 当 丰 >0 时 ， 
Y; =¥,+0A/(t,,¥/)+ (1-0) AS (tY) (3-4) 
Ya =Y; tet Y) AW, +h(t,.%7 AN, (3-5) 
因为 允 是 G3-4) 的 解 , MAY EZ FEIN, 增 量 AN = NG) —N (ty) ， 服 从 参数 为 人 A 的 
HAA: BE AW, EW (ta) -W (t) DAEAR BIRA N (0, A) 的 高 斯 随机 变量 。 
Lice RAIN SEUNG. FR GESTED. Ht efit) th EX 
Y(t)=¥,+0(t-4,) f(t% )+(1-0)(t-4) S (tY). 
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+g (to (P, -W, )+ h(a )(N, =N, ) (3-6) 
现在 将 (3-4) 和 (3-5) 逐 次 代入 (3-6) 可 得 : 

Y(t)=%.+O(t-t,) £(4.% )+(1-0)(t-4) f (hk) 
+g(t Y; )(W,- W, )+A(t Y; )(N,- -N, JHA (ti ¥ a) 
+(1-O) AF (tefja) + 8 (ti Yia) AW, + A(t Ya ANG 


=Y +0(t-1)f (tY; )+(1-0)(t-t)f (tY) 
telak (=W, Jalak (NN Jo af (527) 
1-05a 1)+ Šeli Y; ) AW, +$ hli Y; )AN, (3-7) 
接着 定义 3 个 示 性 函数 
Z (=$ Ttun (Yet O 
Za(t)= YM tan C)+ l) 


N-1 


$=) tpsa) (th 


其 中 加 (D)= i nan 于 是 式 (3-7) 可 改写 为 
Y(t)=Y,+ fo -0) f (5,Z,(s))ds+ fo f (5,Z,(s))ds 


+ [6z (s) dW, + fats, Z, (s)AN, (3-8) 


3.1 近似 解 的 估计 


为 了 证 明 全 局 误差 的 收敛 速率 ， 本 节 将 对 SS6 算法 的 近似 解 作 估 计 ， 引 入 如 下 引 
理 。 


引 理 3.1 设 f :[0,T]xR > R HEBR 3.1 和 假设 3.2, VAY, PERE SSO 算法 得 


到 的 近似 解 ， 则 对 于 任意 的 0<6<1, 4 asmin(l a, 存在 非 负 常数 Ci,C; ， 使 得 


E(k 1%) sce ((% fl} +c 
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证 明 : 由 (3-4) 和 基本 不 等 式 2ab < a +b’, PAN 

lel =i arlag) (0) Gn) 
=P on (at; f -Ar (eX +2007 (1,047 )% 
+2(1-0)Af (t Y,)Y, +20(1-0)A° f (t Y.) f (to Y) 
<P ar (g S -Dar er an, Hasr 
+(1-A) Al, +(1-) | ¢(1,.%)f +00- (4..%)] 
+0(1-0)A?|F(1,.¥%7)f 
<(+A Jf (a)y lr ) + (010)AG+ A (4%) - 
由 假设 3.2 可 知 

lee saraf + (a+ 6a"), (1+ [ef +00) (+4) Ka (1+ 1%) 
<(1+A)(1+(1-0) AK,)|¥,P +20AK, |Y}; +A(1+A)K,. 


HF0<@<1, asmin{ at 所 以 2K,6A <2KsA < 故 移 项 整理 后 得 


2 
2 _(1+A)(1+(1-@)AK,), ，。 A(1+A)K, 
eh! ae Mii) Pe lie a 
1-2K,0A 1-2K,0A 
<2(1+A)(1+(1-6) AK, )|¥, | +2A(1+A) K, 
<2(1+(2K, +1))|%[ +4; 


P| 


两 边 取 条 件 期 望 ， 且 令 C, =2(1+(2K,+1)),C, =4K,， 可 得 引 理 的 结论 。 
引 理 3.2 若 三 :[07]xR 一 R，8g:[07]xR 一 人 和 天 [07]xR 一 人 都 满足 假设 3.1 和 
假设 3.2, YAY 都 是 依据 S$S9 算法 得 到 的 近似 解 ， 则 对 于 任意 0<9<1,， 当 


A smin} fet 存在 非 负 常数 C,,C,， arr ERIO sc. 


=z 
4K, 


WH: 对 (3-6) 两 边 取 平方 后 ， 再 取 关 于 . 务 的 条 件 期 望 ， 得 到 
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E (Yul 156) == (ul L%) += (loos (ox) | a ir |%) +48 (l(a 15) 
+(4A+27A*)E (人 k) yi [ia] +2008 ( (7 (45% A J5) 
+2(1-0)AE((f (1%) 4 )|5)+244E((x,h(r Y; Ja 
+2(1-0)WE((F (tY: ) S (t t,))| 5) +24 E(( F (tre h(t A) 


+24(1-0)a*E(( ny altt A) 
利用 基本 不 等 式 2ab < a? +b, W 
s(t 196) = (Wel 15%) +# (lar (er ls)+ (0-0) er) 
+ AE etri fla) aen) aeg J [5] eae (el ts)+ oa( fia) 
+(1-9) AB (WF Z) + 0-0) AE (| (2) |%) + Az (EFIA) 
if ABE (nl, Y ji [ys] +(1- ojan E|| slar | Is] 
+(1-0) 7 (|F(4,.% )f 156) + aese || serf x) ae 人 (下 I] 
+A(1-0) ME ([F (4% Jf A) +416) ale J LF] 
<(1+4+ 4)E(nP ls)+ar(le(s,n I) 
+((1-0)A+(1-0)a? +(e) 4") EI EFIA) 
+(04+60A? +404? )z| Fey f |] 
+(24A+ 7A? + Aa?) (n(n, yf 四 
由 假设 3.2 可 知 
E( Rul’ |%)<(1+4+ 24+ ((1-0)A+(1-0)A? + 4(1-0) A”) K, E(P ls) 
+(OA+0N? +A0A? +A4+24A4 A707 + ad?) K,E( Kf 5) 
+(2A+ A? +2AA+ AA? +244?) K, 
利用 引 理 3.1 的 结论 ， 可 知 
E (Mal |%)<(+C,a)z(|%[ A) +A (3-9) 
这 里 
C; =(3+4\+7)K, (G, +1), 
Cy = (1+ A+(24+A)K,)+(3 +4441") KC, 
对 (3-9) 式 逐次 迭代 ， 得 到 
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E (al |%)s CAD (+CA) +(14+C,A)" EY} 
i=0 


k+l 


< ((1+C,a)"=1)+(14C,4) E\y,P 
6 


N 
RF es sw A= Talis | <ecz ,所 以 
E((Yal'|%)< = (SEY > (SEY ape 
k+l 0 "O N N 0 


N N 
z% Ca -1 (+ EII 
Č N N 


SC, 


č 
EHC, “eh 


e —1)+ e" E 所 | 
AARSEL ETEN ERa |) SC, BEC, = GG +). 
MEX Y (t)—Z, (t) FY (t)- Z, (t) HHA REET o 
引 理 3.3 设 f:[0,7T]jxR 一 R，g:[0,TjxR 一 了 R 和 h:[0,7]xR 一 都 满足 假设 3.1 和 


假设 3.2， 芭 和 产 都 是 依据 SS6 算法 得 到 的 近似 解 ， 则 对 于 任意 的 0<6<1,， 当 


A < min iea 
4K, 


| 存在 非 负 常数 C,,Cs。 ， 使 得 
E(|7()-2(0 1%) sd, 
E(\¥()-2,(0)f lE) Ga 
证 明 : 当 te[h,tn) 时 ， 由 (3-6) 式 可 知 
Y (t)-Z,(t)=0(t-4,) F(t.% )+(1-O)(t-4) f (tY) 
+ g(t (W, =W, )+h (tole )(M =N, ), 


利用 假设 3.2， 基 本 不 等 式 (a+b+c+d) <40 +4b?+4c?+4d? 可 得 
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E(P- 1%)<40°aK, (1 + E (beef |) +4(1-0) °K, (1 +E(l i 5) 
+4AK, (1 + z(e] |%))+4(4a+ 2202) x, (1 +2(|Kf | L5) 
<4(1-0) AK, E (|¥ y |%)+40%a°K, 
+4(1~0) A?K, +4AK, +4(AA+ XA? )K, 
+ (40°A?K, + 4AK, +4(AA+2A)K,) E(x 5) 

利用 引 理 3.2 可 得 
E(lv ()-Z (0 1%) sca 
HPC, = 4K,C, +4(2+2+A7)K, (C, +1)+4K,. 
HF re[t,,t.,,)B%, RE Z (t)z (t) 的 定义 及 (3-4) 式 ， 可 知 
Z,(t)-Z,(t)= OS (4% )+(1-O) AS (4..%) 
同样 利用 假设 3.2， 基 本 不 等 式 (a+b) <20 +20 WB 
E(|z.(0)-2,(0f 5%)=E(lev (1.27) (0) (a l5) 
<26°A°K, ( +2 (bef Lx) }+20 -6) AK, ( +E(| f L5) 
利用 引 理 3.2 可 得 

E(|2_(¢)-2,(#)/ }%)s Goa 

其 中 C, =2K,(2+C,+C,). FË 

E(r()-20F |) $22 (7 -zoF 1%) +28 (2. -2 (Of LS). 

<C,A 
3X HC, =2C, +2C,. 
3.2 ”收敛 速率 
本 节 将 给 出 SS6 算法 的 收敛 速率 。 


定理 3.1 在 假设 3.1 和 假设 3.2 下 ， Hasmin] fon 存在 非 负 实数 C ， 使 得 方 


-e 
4K, 


#2 G-1) 的 解 X(t) 和 方程 (3-8) 给 出 的 连续 时 间 近 似 解 Y(t) 满足 不 等 式 


华东 理工 大 学 硕士 学 位 论文 第 31 页 
(solr -rof lz ]<ca 
证 明 : 由 方程 (3-1) 和 方程 (3-8) 可 知 
¥()-¥(=(-8) [x e) -7 (8.2 (0) + ol (sx0)-/(5,2())s 
+ a(x(0)-s(52 CO) + fsx (0) -a(5, 2 (yn, 


两 边 平方 后 ， 先 取 上 确 界 再 取 条 件 期 望 ， 可 得 


JrGxO)-7G 


s(a -rof 2) <40-07 z| a 


EU A 
anla [(e(s.X(s))- 5 | 


| (3-10) 


vat anlata 


接 下 来 ， 依 次 对 (3-10) 式 右边 的 四 项 分 别 做 估计 。 
对 于 te [4,tin)， 根 据 Holder 不 等 式 ， 假设 3.1 可 知 


sre n pA 


4(1-0) E 
Osrsi 


t 


<4(1-6) TE 


Co ae 


0 


tied s- ss) 


i 


+8(1-0) TK? zl X(s)-Z,(s)| ds 


<8(1-0) TK? SN (s—t,)ds+ bos 


5 


+16 (1-0) TK; zf fx (s)-¥(s)[ ds 
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由 于 六 ke- ida Je- t, jds ÉM AINA = 一 TA. 再 利用 引 理 3.3 可 得 


: 


<167K2 [E| supl X (r)-¥ (r |Z ds+4T?°K? (4C, +1)A (3-11) 
0srs 
0 PSs 


+16(1-@) TK?E af Ire) AOLA 


2 


4(1-0) E 
O<rst 


Irex- 


类 似 地 估计 (3-10) 右 边 的 第 二 项 可 得 
JexoO-rGzag 网 加 


40 E | sup 
Osrst 


<16TK? Je(sup lx (r)-r (rf L fs +472K (4G, +1)A (3-12) 
0 0<rss 
对 (G3-10) 右 边 的 第 三 项 先 利用 B-D-G 不 等 式 得 到 


letxe Eza] i] 


<isz[ he (s. X (s))- (5.2, (fal) 
再 利用 假设 3.1 和 引 理 3.3 可 得 


JleCxO)-s(zC)jmr 
SGK jz (sup |x (r)-Y(r [% }as+ 67K? (4C, +1)A (3-13) 


| 


4E E 
O<rsi 


ee 
Osrst 


类 似 地 估计 (3-10) 右 边 的 第 四 项 ， 有 


4E | sup 
OSrst} 9 


[(A(s,X(s))-2(5,2, (s))}aN, 


< 中 |(a(s, X (s))-R(5, Z,(s)) WÑ, | 
se al f(a(s,X (s))-h | a 
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di 
中 


+TK?A(4Cs +1)(324 +84°T) (3-14) 


将 (3-11)-(3-14) 代 入 (3-10)， 整 理 得 
(suplx -rE IE] 


<| 167K? +16TK? +64K? +4K? (324+84°T)] ja( sup|X(r)-7(r)| |% je 
0 O<r<s 


t 


<(324+82°7) ) f(a (s,X(s))-A(5.Z, (s))f 


<4K? (324+84°T) 四 可 x (r)-¥ (rf 
0 0srss 


+41°K? (4C, +1)A+47°K? (4C, +1)A 
+16TK? (4C, +1) A +TK? (4C, +1)A(322+82°7) 


运用 Growall 不 等 式 ， 整 理 后 得 
#( sup|x(r)-¥ (rf [%)sca ; 


其 中 


C= gam +64K? +4K} (32448277) |r 
x {47°K? (4C, + 1) + 47 K? (4C, + 1)+167K; (4C, + 1) + TK? (4C, + 1)(324 + 34°7)) ; 
下 面 用 一 个 例子 来 比较 欧 拉 算法 和 SSO BY. WT =4,X, =10,A=20,0=0.9.47 
程 (3-1) 的 系数 分 别 为 了 (t,x) =2x,2(t,x)= 了 (ea = -2x 时 ， 方 程 (3-1) 的 解析 解 为 


L$. 


xof 中 oprim]. 
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图 3.1 精确 解 、 欧 拉 算法 数值 解 和 裂 步 法 数值 解 


Fig. 3.1 Exact Solution, Solution of the Euler-Maruyama scheme and of the Split-step 9 scheme 


图 3.2 离散 解 与 精确 解 之 间 取 对 数 后 的 误差 

Fig. 3.2 Point mean square errors of the Euler scheme and of the Split-step O scheme in a logarithmic scale 

图 3.1 给 出 A=0.01 的 Split-step@ 算法 的 近似 解 和 Euler-Maruyama 算法 近似 解 及 精 
确 解 ， 其 中 近似 解 是 100 次 计算 后 的 均值 。 从 图 3.1 可 以 看 出 Split-step9 算法 的 数值 更 
近似 于 精确 解 。 图 3.2 显示 的 是 Split-step9 算法 的 近似 解 和 Euler-Maruyama 算法 近似 解 
与 精确 解 之 间 的 误差 ， 其 中 误差 是 100 次 计算 后 ， 近 似 解 和 精确 解 差 的 平方 的 均值 。 由 
图 3.2 可 知 , Split-step 0 算法 的 近似 解 与 精确 解 之 间 的 误差 更 小 。 此 例 说 明 Split-step9 算 
法 要 优 于 Euler-Maruyama 算法 。 
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第 4 章 ” 带 跳 的 随机 延迟 微分 方程 2 法 近似 解 的 收敛 速率 
本 章 在 上 章 的 基础 上 ， 在 方程 中 引入 时 滞 项 ， 考 虑 如 下 带 跳 的 随机 延迟 微分 方程 : 
aX (t) =f (t, X (t), X (t-r))dt+g(t,X (t), X (t-1)) aw, 
+h(t,X(r),X((t-2) ))aN, rs[o7] (4-1) 
X (t)=¢(t) t e[-z,0] 
其 中 到 为 .9 适应 的 一 维 布朗 运动 ，NN, 为 参数 为 4 的 一 维 泊 松 过 程 ， 且 到 与 ,相互 独 
立 。X( 太 )=lippXX(s)。 对 于 所 有 的 1>0,s e[-r,0]，g(t)eC[-r,0]，p(s) 与 所 和 WN, 相 
互 独立 , E supa lp(s) <0 GD) 为 属于 C([0,7];R) 的 随机 变量 , 初 值 为 (0) = 9(0) » 


7 为 一 固定 的 正 时 滞 . 假 设 方程 的 系数 f,g,h 为 其 所 有 变 元 的 连续 函数 。 
为 了 保证 方程 (4-1) 存 在 唯一 解 ， 还 需要 引入 如 下 假设 
假设 4.1 系数 f,g,h 关 于 t 满 足 Holder iE, KFA x, y WA Lipschitz 条 件 ， 即 


对 于 任意 的 0<t <7,0<4,5TWURx,x,,y,y,€R, FEEBRK, > 0， 使 得 


UEBA CESA M CERA CEA KA ETA a CEASA) 
1 
SK, (k -bl +|x, -x+y -»,| (4-2) 


假设 42 系数 f,g,h 均 满足 线性 增长 条 件 ， 即 对 任意 的 te[0,7] 以 及 x,yeR， 存 
在 正常 数 玉 , > 0， 使 得 
[fx yeer vex <.K,0 +x) +f) (4-3) 
假设 4.3 ( 初 值 连续 性 ) 对 于 任意 的 s,te[-r,0]， 存 在 正常 数 B8， 使 得 
Elp(s)—e(of <4|s—t 


由 文献 59 可 知 ， 方 程 (4-D) 存 在 唯一 解 。 
现在 ， 本 章 将 对 带 跳 的 非 线性 随机 微分 方程 (4-1) 建 立 SS68 BIKE. MAAK, WE 


T=NA, T=mA, 其 中 N,m 均 为 正 整 数 .节点 t =kA, k=—m,...,0,1,..,N . 当 大 <0 时 ， 


¥,=y(kA); TH4k >On, 
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yy =F HOME (Phim) +(1-O)AF (Yo Yen) (4-4) 
Ka = H+ 8 (tole Yn) AW, +A(t e JAN, (4-5) 
其 中 增 量 AN,= NG。)- NG:) ， 服 从 参数 为 4A 的 泊 松 分 布 ; HEE AW, EW a) -W A 
相互 独立 的 且 服 从 N(0,A) 的 高 斯 随机 变量 。 

对 于 上 述 离散 化 近似 解 二 ， 构 造 时 间 连 续 的 近似 解 。 对 任意 的 te [6,t,,)， 定 义 
¥(t)=¥,+0(t-t)F (to Yi Yim )+ (1-9) (0-4) f(t, FY) 

+8 (te Fes Yen )(W,-W, )+h(to Yen (N,N, ) (4-6) 

与 第 三 章 类 似 ， NE RONG -6)， 整 理 可 得 
Y(t)=F +O(t-t) f (4.97 ¥en) +(1-0) (t-t) f (too Yrm) 

a ewe ape 


+(1-0) 5 Af (t, kn) + s(t Wn) AM +$ h(t Yn JAN, (4-7) 
接着 定义 3 个 示 性 函数 


N-1 
Z (t)= 2 Y Xunn) (C)+ Vv Leary (t): 
N-1 
Z,(t)= ZY Ay sa) (C)+ Fv Ze (t), 
- N-l 
a 2 1X dia) (7). 


其 中 入 (‘) ARERR, M4rer HM, y(t) Bil g (1)=0。 于 是 式 (4-7D) 等 价 于 


Y(t)=%+ Ja- f (5,Z,(s),Z,(s—mA))ds + fo f (5,2, (s). Z, (s—mA))ds 


+ fel Z,(s),Z,(s—mA))dW, + fr (5,Z, (s),Z,(s- mA))dN, (4-8) 


4.1 近似 解 的 估计 
为 了 证 明 全 局 误差 的 收 敏 速率 ， 本 节 将 给 出 一 些 引 理 。 首 先 对 到 (信人 | 现 ) 和 
B(ly j z) 作 估 计 。 
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引 理 4.1 设 f:[0,T]xR 一 满足 假设 4.1 和 假设 4.2, 和 瑟 都 是 依据 SS8 算法 得 


到 的 近似 解 ， 则 对 于 任意 的 0<2<1， 存 在 非 负 常数 Ci,C。， shcmin| hat 有 


8K, 
max E ( A 


Osrsk 


|%) sc, max E (|Y, Í 由 jc 


Osrsk 


证 明 : 首先 给 定 某 个 整数 ke[0,N]， 对 于 任意 的 非 负 整 数 re[0,k] ， 运 用 引 理 3.1 


同样 的 证 明 方 法 ， 先 对 (4-4) 式 两 边 取 平 方 ， 再 利用 基本 不 等 式 2ab <a’ +b? ARK 4.2, 
得 到 


Iyf <(1+4) a 


<(1+A)(1+(1-4)AK,)|¥,[° ol P Ear 
+(1+4)0AK, tai 上 ')+a(1+4)K, 
上 式 的 两 边 取 条 件 期 望 得 到 
人 | 醒 jsG+A)G+G-e)Ak)z(z |%)+(1+4) (1-6) AK, E (Ent |-%) 


y, 


* + (0+0?) K, (1+ + E) 


F-n 


+(1-@)A(1+A)K, (1+|¥,) + 


+(1+A) OAK, (z( 


rila) z(e: f |%))+a(+ 4) x, 


< (1+-A)(1+(1-0) AK, max £ ([¥/']%) +(1+4)(1-0) AK, max E (lf 15) 


Osrsk 


+(1+A) OAK, (max E(I |%)+max (le) 四 外 A(1+4)K, 


<(1+A)(1+(1-0)AK,)max E(P |%) 


+(1+A)(1-0)AK, (max B(\¥._.['|%}+ Elet) 


msrsk 


+(1+A) OAK; (max E (rf |) +max (lr 


Osrsk r 


站 二) Elof }+ A(1+A)K; 
<(1+A)(1+2(1-6) AK, )max E(P I5) 
+2(1+A) OAK, max E(|¥"['|5%)}+A(1+4)K, (1+ Elol ) 
注意 到 上 式 的 右边 与 x 无关， 于 是 有 
max EE |%)<(1+A)(1+20-0)AK, )max & (|r [ |) 


Osrsk 


+40AK, max (rT 55) +a +a) K (1+ Eol) 


0<0<1 和 A<min 人 二 | ; PLLAKA <4K, AS, 故 移 项 整理 后 得 


2 
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2 (1+A)(1+2(1—@) AK, ) AQG+A)K 2 
mag Eff 56) a Baas (al I) ae (1+ Ee) 


<C, max E(|¥,['|%)+C, 
RPC =4(1+2K,),C, = 4K, (1+ Efe) 
引 理 4.2 设 f:[0,7]xR>R, g:[0,7]xR>RAA:[0,7|xR>R BREE 4.1 和 


假设 4.2， ALY, ABE KE SSO 算法 得 到 的 近似 解 ， 则 对 于 任意 的 0 < @ <1, 存 在 非 负 党 


数 C,,C,， MAMA smn} that 有 
8K, 


max E ( 下 
Osrsk 


gea max E |E |%)s C,. 


Osrsk 


证 明 : 当 0<k<N 时 ， 由 (4-6) 式 可 知 


区 中 = 区 ,到 中 +- ion) +e (iYi) Amf 


2 
+|n( tJ, 


Ye »)AN,| +2¥A(t5%; Yim) AN, 

gover 六 )¥, +2(1- der Yam) Y +2,8 (tY Yin) AW, 
+2(1-0) 6A? Ha ERE A We RE 
+20M (bo Yg Yin )h(t Xin) AN, +2(1-O) AF (tsk, Ym) 8 (iYi ) AW, 
+2(1-0)Af (ty Yeo Yrm) h(t Y; Yim) AN, +28 (tof; Mein AM A(t Y; Yin) AN 


AH AW,, AN, 和 ,多 Kit, BB) EAW, Z)=0,E(AN, | Z)=4. ERE 
32 相 
类 似 ， 对 上 式 两 边 取 条 件 期 望 ， 运 用 基本 不 等 式 和 假设 4.2， 得 到 

E( Meal ta) aaae E) eleta) Elta) 
+((1-8)A+(1-0)4? +4(1-0)A7)K :+ 人 |) +z(| Mi z) 
+(0A+0A? + 20M?) K, (elf La) e(l] Ea} 


+(24A+A7A? + AA’) K, ab Aa = 加 
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其 中 CG =(1+2)+(2+A)K,,C, =(34+44447)K,,C, =(2+4)K,. 
不 等 式 的 两 边 取 期 望 。 与 引 理 4.1 相 类 似 ， 给 定 某 个 整数 ke[0,N]， 对 于 任意 的 整 
数 ye [0, k] ， 有 
(Ivo }%) = (1GA) may (eh a) + Ca (ma (P|) + e(a) 
+ČAmax E(|Y n|) +24 
<(1+ČA+Č,A)max E(|Y, f (5) 
+CA (mas Er La) + mas E(f) ENa) 
+GA (ma slr |%)+ Elo) +E 
<(1+ČA+Č,A) max (|v |.) +2C,Amaxé (| L5) 
+(G +2 )a(1+ Eloi) 
利用 引 理 4.1， 可 得 
Yal |%)S(1+C,A)max E(P A) +CA (4-9) 


a 


其 中 CG =C,+2C,C, +C,,C, =(Č, +C,)(1+ Elf )+2C,¢,. 


注意 到 上 式 的 右边 与 + 无 关 ， 且 对 于 任意 的 x e[0,4] 成 立 ， 故 有 


max E(X a} |5) < (1+C,A)max E(f 1%) +C.4, (4-10) 


Osrsk Osrsk 
对 (4-10) 式 进行 逐次 迭代 ， 得 出 


k+l 


k ' 
max E(f af 1%) < CAD (14 C,A)'+(14C,A) E|Y,P 


<S4((1+C,a)"-1)+(1+6,4)" Blof? 
用 引 理 3.2 相同 的 方法 ， 得 到 
max E(|y,.,[ |%) <C, 
O<rsk 


C, 


EFC, a all —1)+e%E|y,P. 
5 


2 
* 
7 


再 利用 引 理 4.1， 便 可 得 出 max E( 


男 ]<C。, 其 中 G0 


引 理 4.3 设 f:[0,7]xR 一 R，g:[0,T]jxR 一 R 和 4h:[0,T]xR 一 都 满足 假设 4.1 和 


假设 4.2, Y, AY, 都 是 依据 SS6 算法 得 到 的 近似 解 ， 则 对 于 任意 的 0<69<1,， 存 在 非 负 
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实数 C,,C,， was sin} gp | 有 
8K, 

E(|r()-2(0 |) SGA, 
E(|7()-2,(0 4) sd. 

HERA: 当 te[ 太 ,大 时 ， 由 (4-6) 式 可 知 
¥(t)-Z,(t)=0(t-4) f (to Yi Men )+(1-9) (0-4) F(4.% ) 

+g (tot; en) (WW +A (40% Yen (N, -N, ) 

5 5 Æ 33 的 证 明 方 法 相似 。 利 用 假设 42 和 基本 不 等 式 

(atb+e+d) <4a +452+4c2+4d2， 可 得 


ern) saree) >a) 


1a) 


+4(1-6) A°K, (1+ r(x | |7)+E (is. 
+4AK, (1 + E(\Wef | %} +2 (lien) 1% )) 

+4(AA+22A7)K, (1 +E ( [He z) E(t] L) 
<4A°K, (e(l [ L% J+ max E | Kal |% J+ E lol’ }+4a°K, 


msksN 


+4A°K, +4AK, +4(4A+27A?)K, 
"| + max E(f a) zlol ) 


msk<N 


+ (447K, + 4AK, +4(4d+ 220°) K,)( (lr; 
结合 引 理 4.2 的 结论 可 得 

E(\¥()-Z,(0f |%) SGA. 

其 中 

C, = 4K, (2C, + Blof’ )+ (4K, +4K, +4(2+2) K,)(2C, + Ellol)+4K, +4K, +4K, +4(4+4?)K,. 
再 根据 (4-4) 式 ， 可 知 

Z, (t)-Z, (t) = OAS (tY; Yin )+(1-0)Af 到 ,及 


同样 利用 假设 4.2， 基 本 不 等 式 (a+b) < 2a’ +2b 可 得 


Sl ict ha — _ MIR 
E(|2, (t)-Z, (Of L%)= E (ayy ( tales Ye») +(1- A) Af (tY, ke | l5) 
<208E(|/( 2 "|%)+20-0) are (|s( 本 /本 | |] 
<26°A°K, (1+ E( (re %) +2 (Heol 1% ) 
V) 
<2A°K, ( + E| 5%} + max E (sl IF)+E le) 
+20°K, (142 (|v, 1%)+ max E (Enf 15) Elo) 
运用 引 理 4.2 可 得 
E(|Z, O-A E) Gud, 


其 中 Ci = 2K, (1+2C, +E |g )+2K, (1+2C,+B|of). 于 是 


+2(1-0)' A?K, (+ 人 fla) 


AORA AE ORAO L%) +22 (|Z ORAA 
<$2C,A+2C,,A 
a ae 


ZEG = 20, +20,,. 


4.2 ”收敛 速率 
BRA SSO 算法 的 收敛 性 。 
定理 4.1 在 假设 4.1， 假 设 4.2 和 假设 4.3 下 ， 存 在 一 个 非 负 实数 C ， 使 得 方程 (4-1) 


的 解 X(t) 和 方程 (4-8) 给 出 的 连续 时 间 近 似 解 了 () 满 足 不 等 式 
(sp x) OPS }s CA. 

证 明 : 证 明 方 法 与 定理 3.1 相似 。 首 先 由 方程 (4-1) 和 方程 (4-8) 可 知 
x()-7()=0-0) (sx 0), XG-mA) -5,200),2(s -mA) ys 


t 


+0 (A( f(s,X(s),X(s—mA))- f(s, Z,(s),Z,(s—mA)) as 
e] g(s, X (s), X (s—mA))-2(5,Z,(s),Z, (s-mA)) WwW, 
+ Q s,X (s), X (s-mA))-h(5,Z,(s),Z, (s—mA))\dN, 
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两 边 同 时 平方 后 先 取 上 确 界 ， 再 取 条 件 期 望 ， 可 得 
s) 


z( sup|x -rof LE 


y 2 


|(7 (s.X (s),X(s-mad))- f (5,2 (s),Z, (s—maA)) ds 


4(1-0) {oe 
Osrst 


y 2 


f(s (5-X (s), X (s-mA))- f (5, Z, (s),Z, (s—mA)) Jas 


0 


+40°E | sup 
Osrst 


+48 fle (s,X (8), X (s-mA))- g(5,Z, (s),.Z,(s—mA)) am, x 
+46 flo(s,x(s) x (e-ma)) -A(z (s),Z,(s—mA)) dN, s 


(4-11) 
接 下 来 ， 依 次 对 (4-11) 式 右边 的 四 项 分 别 做 估计 。 


ME elt sta)» WH Holder 不 等 式 ， 假 设 4.1 可 知 


4(1-0) E i Ff (s,X (s), X (s-mA))- f (5, Z, (s),Z, (s—mA)))ds 


2 
z) 


中 


<4(1-0) T [e(l xl), (6-ma))~F(,24(6).2(6-ma) 


1 1 i 
<12K}T | fx(s)-Z,(s)f al | 
0 0 


9 


由 于 flo—sfas= Si [let )do+ f(s- J< ETA RBI 43 和 假设 4.3 可 得 


2 
s 


+ pr [pe om -Z,(s—mA) ds 


r 


fl (sx (8), X (s-ma))- F (82 (8), (s-mA))jas 


0 
s) 


4(1-0) | sup 
Osrst 


t 
< scrariakzl fx(s)-2,(s)f as 
0 


narr 了 |X (u)-Z,(u)f du 
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<6K272A+ prers[ fx (s)-Z,(8)f as s) 
0 


4 


< srrat pirs[ Ix()-2 (s)f as 
0 


+ wax T (u)—Z, (u)f du 


| |X (u)-Z,(u)} du 


s) 


s 


<exir'as asxar( fe sup [x(r)-¥(r)f 
0 Osrss 


+12K7T BTA 


t 
<6K272A+ saxire| Nx(s)-7(s)| as 
0 


+12K?T Br +47 re fir(s) -Z, (sj ds 
0 


z) as s10K pa +48K°TC,A 


(4-12) 
类 似 地 估计 (4-11) 的 右边 第 二 项 可 得 
ce sp (1 (s.X(s),X(s—ma))- 7 (8, Z,(s),Z, (s—mA))}ds 加 
srst 0 
<6K2T2?A+ ‘sper fel aw |x(r)-¥(r)f 5a) 
+12K?T BrA+ 48K2TC,A | (4-13) 
对 于 (4-11) 的 右边 第 三 项 利用 B-D-G 不 等 式 得 出 
可 到 [le (s, X(s),X(s- mA))- g (5,2, (s),Z,(s -mA)) WW, g 
Osrst 0 


<16 [a(|e(s,X(5),X(s—mA))-g(5,2, (s),Z(s—ma))} L% xs 
再 利用 假设 4.1 和 引 理 4.3 可 得 
| 


Osrst 


r 


[(e(s.x(s).x (s—mA))-g(5,Z,(s),Z,(s- mA))\aW, 


0 
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<24TK7A+192K? fa sup [x (r)-Y(r)f B 
0 Osrss 


+ 48K? BrA+192K2C,A (4-14) 
类 似 地 可 得 

4E (æ J (s), X (s-mA))-h(5,Z, (s),Z, (s-a) Is 
i E fals (s). (s—ma)) -a(s Z,(s),Z,(s—mA)) aN, 网 
ol 下 f(r(s.x (s),X (s—mA))—h(5,Z, (s),Z, Gn) | 


<(324+84°T) fe (os. (6). X (6—ma)) -nz (s),Z, (s—mA)) | fas 


3 


< 5 TKI A(324 +81°T) +12K? (324+82°7) J E ( sup [x (r)-¥(r)f Ls es 


+12K7C ,A(324 +8A°T) + 38K; A(324 +8277) (4-15) 
将 (4-12)-(4-15) 代 入 (4-11)， 整 理 得 


a( sup [x (r)-¥ (rf lz) SCG fe | sup [x (r)-¥(r) Ls 十 CA 


其 中 Ci =96TK? +192K? +12K? (32448277), 
3 
2 

+48K?TC, +48K27C +192K?Cs +12K2Cs(324+8227)- 
运用 Growall 不 等 式 ， 整 理 后 得 


Č, =1272K2+24TK2 +—TK? (324+ 84°T)+24K/T Br + 48K; Br + 3K; Br(324+82°7 ) 


#{ sup|x(r)-¥ (rf GE ČA, 


其 中 C=Ce5. 


由 于 sup [xX (r)- ¥(r)f = sup |x (r)- ¥(r)f v sup [x(r)- Y(r)f ， 运 用 假设 4.3 得 
一 TSrS1 -TSr<0 Osrst 
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B{ sup x (7) -YF 1%) E(x- r E)E, su Ix(r)-7(r) | 


0 人 


<CA+E 


<CA, 
其 中 C=C+B。 
例 4.1 讨论 下 列 带 跳 的 随机 延迟 微分 方程 
dX (t)=|aX (t)+bX (t-1)]dt+cX (t-1)aW (t)+ DX (t)dN (t),t20 (4-16) 
. X (t)=1+1,t e [-1,0] 
对 式 (4-16) 的 系数 取 a= -9,b =7,c =1,D =0.2,A=1,0=0.9. 精确 解 是 根据 [38] 的 方法 得 到 


的 , 即 取 步 长 1/1024 的 Euler-Maruyama 算法 的 近似 解 作 为 精确 解 , 图 4.1 给 出 Split-step 
9 算法 的 近似 解 和 Euler-Maruyama 算法 近似 解 ， 从 图 4.1 可 以 看 出 Splitrstep 0 算法 的 数 
值 更 近似 于 精确 解 . 图 4.2 显示 的 是 Split-step9 算法 的 近似 解 和 Euler-Maruyama 算法 近 
似 解 与 精确 解 之 间 的 误差 ， 由 图 4.2 可 知 ，Split-step 算法 的 近似 解 与 精确 解 之 间 的 误 
差 更 小 。 此 例 说 明 Split-step9 算法 要 优 于 Euler-Maruyama 算法 。 


| 
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图 4.1A =0.01 BY Split-step 算法 、Euler-Maruyama 算法 和 精确 解 模 拟 结果 
Fig. 4.1 Exact Solution, Solution of the Euler-Maruyama scheme and of the Split-step O scheme with 
A=0.01 
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图 4.2 离散 解 与 精确 解 之 间 的 取 对 数 后 的 误差 
Fig. 4.2 Point mean square errors of the Euler scheme and of the Split-step @ scheme in a logarithmic scale 
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第 5 章 带 跳 的 随机 微分 方程 裂 步 法 一 阶 近似 解 的 收敛 速率 


本 章 考虑 带 跳 的 随机 微分 方程 如 下 
aX (t)= f (t, X (t))dt + g(t, X (t))aw, +h(t,X (F ))aN, 
X(0)=9(0) 
RAW AZ 适应 的 一 维 布朗 运动 ， NABMA ANH, BW, 5 NAR 


(5-1) 


立 。 X(C)=limX (s). WER >0, 90) SW, M N, 相互 独立 , 且 Elp(O)| < o0; 


X(D) 为 属于 C(0,71R) 的 随机 变量 。 假 设 方程 的 函数 gh ARRATERA, 


关于 第 一 个 变 元 存在 一 阶 连 续 偏 导数 ， 关 于 第 二 个 变 元 存在 二 阶 连 续 偏 导数 。 
现在 给 出 本 章 的 假设 条 件 : 


假设 51 系数 /,g,h 关 于 第 一 个 变 元 具有 一 阶 连续 偏 导数 有 界 ， 关 于 第 二 个 变 元 
的 一 阶 偏 导数 存在 且 连 续 有 界 ， gg', gi 满足 Lipschitz 条 件 。 

假设 52 系数 /,g,h 均 满足 线性 增长 条 件 ， 即 对 所 有 的 0<1<T 以 及 xe R， 存 在 
正常 数 K, > 0 ,使 得 [Fez veeo vjh(i,x) < K+|x| ). 系 数 f,g,h 关 于 第 一 个 变 元 


的 一 阶 连续 偏 导数 满足 线性 增长 条 件 ，g? /满足 线性 增长 条 件 。 
由 文献 的 可 知 ， 在 假设 5.1 和 假设 5.2 下 方程 (5-1) 存 在 唯一 解 . 
在 第 三 章 中 我 们 给 出 了 方程 (5-0) 的 一 个 收 敏 速 率 为 的 SS9 算 法， 与 第 三 章 不 同 ， 


本 章 我 们 将 通过 tô- 7ayyor 展 开 式 构造 一 个 收敛 速率 为 1 的 Split-step 算法 。 
5.1 基于 116 一 Taylor 展开 的 Split-step 算法 


首先 给 出 半 拷 116 公式 。 假 设 1(/,X(1) 关 于 1 一 阶 连续 且 关 于 x 二 阶 连续 ， 运 用 半 扶 
Hô 公式 ( 见 文献 [3]) 知 : 


1 1 


L(t, X(t))=7(0,X(0))+ f Lal (s, X (s))ds+ f LI(s,X (s)) aw, + Lilo x(o)an, (5-2) 


0 0 


其 中 
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=1 (t,x)+ f(x) (tx)+58" (sax) (t,x) (5-3) 
LI (s,x) = g(t, yh) = g(t,x)l,’ (t,x) (5-4) 
L_l(s,x)=1(t,x+h(t,x))—1(¢,x) (5-5) 


在 给 出 116 一 Taylor 展开 式 之 前 , 我 们 需要 引入 一 些 定义 和 符号 。 由 于 在 116 — Taylor 
展开 式 中 涉及 关于 时 间 t 的 积分 ， 关 于 维 纳 过 程 成 的 积分 和 关于 泊 松 过 程 的 积分 ， 为 了 
叙述 简便 起 见 ， 本 章 将 分 别 用 0,1,-1 来 表示 上 述 不 同意 义 的 积分 。 由 于 在 推导 过 程 中 ， 
一 般 都 会 将 泊 松 过 程 转化 为 补偿 泊 松 过 程 ， 所 以 我 们 用 一 1 来 表示 补偿 泊 松 过 程 的 积分 。 
又 由 于 116 一 Taylor 展开 式 中 含有 多 重 积分 ， 因 此 我 们 需要 引入 多 维 指标 。 定义 一 个 行 疝 
Bo=(j,ji), BP FE{-L-101}, 1<i<J,JEN, RX MIMBAS BBE 
wa 。 这 个 多 维 指标 的 长 度 .Je {1,2,3,....} 0 AM Ip) =a = 二 2。 为 了 使 这 个 定义 更 
为 完整 , 用 表示 长 度 为 0 的 指标 , BJ, =0, 这 意味 着 不 对 函数 进行 任何 积分 。 另 外 ， 
DAV p, n, Aw, RRS Bia F- 0, 1A. HH 


Po = 95 Mo) = 1, Wo = 0, Pao) = Qro) = L Wao = L Pay = ba-y = 9 Way =1- 
将 所 有 多 维 指 标的 集合 定义 为 
w= {v}Ufa=(A i. d)):4, €{-L-LOM1<i< J.J EN} (5-6) 


我 们 定义 一 个 不 含有 -1 多 维 指标 集合 ， 即 
p={aep:je{-1,0,1},1<i<J,JeN} (5-7) 

非 空 集合 A c A AMEN B(A)=f—A. 

给 定 aefJ,>1， 对 于 及 上 的 元 素 ， 用 一 a，a 一 分 别 表示 删除 多 维 指标 的 第 一 
个 元 素 或 最 后 一 个 元 素 ， 例 如 
= (jo doses) = (jz) (5-8) 
a= (J, jans ji) = (hja) (5-9) 


又 例如 ， 若 a=(1,-1)， 则 一 a= 一 (1,—1)= (—1),a+=(1,-1)-= (1). 
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最 后 , 对 于 户 上 的 任意 两 个 多 维 指标 = (jojo dy) B= (ihi) 定义 算 子 入 
的 运算 为 
BRO =i ETE TS hse) = Gielen lathes 
例如 ,对 于 a=(1,-1) 和 6=(1,0)， FaAP=(1,—1)A(1,0)=(1,—-1,1,0)- 
下 面 给 出 rô- Taylor 展开 式 。 
在 假设 5.1 和 假设 5.2 下 ， 方 程 (5-1) 的 解 X(1t) 可 写 为 如 下 等 价 的 积分 方程 形式 


X (t)=X (0)+ [f(s,X(s))as+ fe(s,X(s)am, + fa(s,X(s))aN, (5-10) 

令 1:[0,T]xR 一 R， 为 一 个 连续 可 微 的 函数 。 当 1(t,x)=x 时 ， 运 用 半 蒜 116 AX, 
得 Lx= 了 ,Lx=g,Lix=h 则 由 (5-2) 式 得 到 的 即 为 (5-10) 式 。 

41=f,l=g,l=h, 应 用 (5-2) 式 得 到 最 简单 的 116 一 Taylor 展开 式 ， 即 
X (t)=X(0)+ | f (0,X (0))+ (it, Ff (r,X(r))ar+ jz f (r,X(r)) aw, + {rs sex (an, Jes 

' [2(0.x(9)+ [te(r.X(r))ar + fhe(rX (r)a, + JrsCxOjav a 

+ roxo) [thlr x (r))ar + [an(r, x (rn))am, + [ar X(r)) an, Jax 

=X (0)+ f (0,X (0)) fds + g (0, X (0)) far, +h(0,X (0)) JaN, +R (5-11) 


余 项 
R= [[l,f(r.X(r)) drds+ fast X (r))dW,ds+ Jfar r,X(r)) dN,ds 


+ 


jms „X (r)) drdW, + J js 人 r,X(r))dW,dW, + J | g(r,X(r))aN,am, 
+ ea (r,X(r))drdN, + J fan(r, X(r)) Wan, + | fral r, X (r)) aN, dN. 


对 上 式 继续 使 用 半 蒜 116 公式 ， 依 此 类 推 下 去 ， 得 到 更 高 阶 的 116 一 Taylor 展开 式 。 
根据 文献 [40]， 对 上 面 的 多 重 随机 积分 ， 可 递归 的 定义 
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X if =Ù (5-12) 
¢ 
[Le[x,at J, >0,j,=0 (5-13) 
P 
¢ 

I, {xf =4 flee [x], dF, J, >0,j,=1 (5-14) 
P 
¢ 
[Lf] aN, J, >0,7,=-1 (5-15) 
Pp 
5 
fle [x], dÑ, J, >0,j,=-i (5-16) 
P 


其 中 en0O<p<s<T，N = WN, 一 At 为 补偿 泊 松 过 程 。 
对 于 ae 朱 ， 定 义 函 数 1:R1xR 一 怀 如 下 : 
x :QC=y (5-17) 
ce ae : ae fi\{v} (5-18) 
这 里 算 子 工 用 式 (5-3)-(5-5) 来 定义 . 
将 (5-12)-(5-16) 和 (5-17)-(5-18) 代 入 (5-11) 可 以 得 出 


XY [z (0,X(0))], + D li (GZ 站 (5-19) 
aeA\{y aef(A 
Kh p={aep: j,€{-1,0,1},1<i<J,0<J<M}, (5-20) 


其 中 M ATERRM. ANRSELER, B(A)=f—-A, PKS-19)H 6—Taylor 展开 式 。 
为 了 建立 一 阶 近 似 解 ， 将 时 间 等 间距 分 割 ， 取 A 为 步 长 ,满足 7=NA ， 其 中 入 均 
为 正 整 数 .节点 =kA ，k=0,1,.…, .根据 文献 是 ` 知 方程 (5-1) 的 精确 解 X(z) 的 


J16 — Taylor 展开 式 为 


xoxo D [Prhe +l) 


acA\{v} \ i=0 


t2, Sa. CX" tlh I (5-21) 


ae f(A) \ i=0 
Hip Ac fi HAE, B(A)=fi-A. 
4M =1 时 ， (5-20)s% A={v}U{(-1),(0),(1)}, TUX Split-step 算法 : 
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Y =y [iy Xe) |” oto tn) (622) 


Ya =E +l (tX o aaler] So 
利用 (5-5) ~ (5-12) 、(5-13)、(5-17)、(5-18) 可 得 


Ty [ia (4.5% i = ls [x P So Ë, 


i 中 = je)a- Af (tY), 
Iliw (4%) = fel Ue Y; )aw, = g (hs )AW,, 


talta ni | = Ta tys¥ JAN, =h(t, Y; )AN,. 
tk 


于 是 (5-22) 和 (5-23) 等 价 于 
Y =Y, +Af (tY, ) (5-24) 
Vn =Y; +g(t,, Y; )AW, +h(t,,Y; )AN,. (5-25) 


BD AREKHERY Split-step 算法 ， 其 收敛 速率 为 Z i 

当 取 (5-20) 式 中 的 M =2 时 ， 
A={y}U{(-1),(0),(1),(0,0),(0,1),(0,—1), (1,1),(0),(,-1),(-1,-),(-0),(-L1)}, 232 
本 章 所 要 讨论 的 基于 116 一 Taylor 展开 的 Split-step 算法 : 


R=%+ DL) (5-26) 


p 


Y,a =Y + 要 LLX e > (1%) |" (5-27) 


wep(A,) 
a=} OED + 8(As)={(41).(4-1),(-11),(-L-D} > 
(5-3)-(5-5)、(5-12)-(5-16) 和 (5-17)-(5-18) 可 得 


Ilol ) |" -jz 人 tY; dW, = g(t% ) AW, 


kaltak) Ji = Ta tY; ) aN, =h(t,,%; )AN,, 


ty 


hala (o¥) = = jel)e (6 awa, = e( Mas (ao) am) -二 


te hy 
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nota lo = Telak Je (oo Jaa, = (XA (on) S Mar} 
talat Tie (gi +n(t,.¥/))-2(a.%¢))an.am, 
(eta -se aman, $ (ra =m) 
Toa [ten (eo) -] Mole +h(t¥))—h(t..% ))aN,an, 
-人 .+ 全 -az SAN.) -Zan } 
这 里 了 (让 代表 第 i 次 泊 松 跳 的 时 刻 。 于 是 (5-26) 和 (5-27) 等 价 于 


Y, =Y, + f (toY,)A, (5-28) 
Ya =Y +g (th AW, +h(ti Y? )AN, + (oY); (tY SCA -A) 
+ g(t h (iF)x To bh +(2 (08 +4(4.¥))-28 (t8 fag LT 
+5 (a(t¥ +h(1,5%))- h(t,,¥;))((ANe) - AN, ). (5-29) 
将 (5-28) 带 入 (5-29) 得 


Ya =Y, + f (to Y.) Atg (tY )AW, +h(i,.¥ JAN, +58 (to). (4,,¥°)((Am) -4) 
+ g(t Fe)h (to Plyg + (g(t +4(4%))-2(4.%))xZay lh 
+5 (n(t¥ +h(t,,¥;))-h(t,.¥;))((aN,)' -AN,) 


=K Jn [the rl"+ > l[i (to y) +5 nuir -o 5630 


acp) 

为 了 定理 证 明 更 为 简便 ， 令 A=A,UA; ={v}U{(0),(1),(-1)}, #EB(A)=A\A, 
Bi B(A)={(141),(1,-1),(-1.1),(-1-1),(0,0),(0,1),(0,-1),(1,0),(-L0)}- E 
B(A)=B(As)VA(As)» &A(A,)={(0,0),(0,1),(0,-1), (1,0),(-10)} « 

注 . Sac A={v}U{(0),(1),(-l} MN, 1 PMSF x, fgh, HER 5.1 可 知 1 满足 
lipschitz 条 件 ; 24 a€6(A,)={(1,1),(L—1),(-L1),(-L—-1)} 时， 有 分 别 等 于 
gg, .gh,,g(t.x+h)—g(t,x),h(t,x+h)—A(t,x), KEBE 514 grh AR ght 
lipschitz 条 件 ， 所 以 ggl gh! 分 别 满足 lipschitz 条 件 ， 随 后 依据 假设 5.2 可 知 ， 当 
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aeAUH(A;) 时 , 1 满足 线性 增长 条 件 ; ae 6(As)={(0,0),(0,1),(0, 一 1),(1,0),( 一 1,0)} 时 ， 


7 


I 分 别 等 于 fl +l +587 Sigl + fer +5 ~ gg" h + fh, 十 二 TURE 


f(tx+h)-f (tx), HF gg! = (gg), —(8) gh! = (gh), -gih RIE 5.1 可 知 
ggi gh 有 和 界 ， 再 根据 假设 5.2 知 1 满足 线性 增长 条 件 。 


5.2 Split-step 算法 的 收敛 速率 
本 节 我 们 将 证 明 Split-step 算法 (5-28) 和 (5-29) 的 近似 解 收敛 于 方程 (5-1) 的 精确 解 ， 
并 给 出 收敛 速率 。 为 了 叙述 简便 起 见 ， 我 们 引入 时 间 连 续 的 近似 解 ， 当 te | 六 ,As) 时 ， 
定义 
¥()= ly [ly (o%)], +o [ho G57)] +o ler), 
+Iey[len (tek, )] + In (tY, ‘I, + Tan [toon (tF; “15 
+ Kaa [ien (tof; K + lo le (oY )], 
= Toy ÜE hoy (th) + oll hoy (40% )+ Ay LT, x (45%; ) 
Hy [i x (oY, Jani, aay (tY, a DY, uy (4 ¥/) 
Hean I, ham (oF) Hy LY, ka (4-27) 
Ke BL) + DL 


K+ frehel Y Jaw,» [fe (1.0!) as (aY, )amam, 


在 在 


+ [fe tY; )h, (tY; WAN, 


iy hy 


+ [s(n salax ))-8 (xt) ava 


+ falt) JaN, mG (to; +h(t, „Y; ) IN, aN, (5-31) 
将 (5-30) 式 逐次 代入 (5-31) 式 ， 可 得 
Y()=%+ > Poa kt" + La [la (te, 


acd, y} 


Ife ab wrt ntr] 


acA(A;) | i=0 


+ 5 1 SA ARA Ji +L lla (4%) Ai (5-32) 
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其 中 A ={v,(O)}, A, ={(1),(—}» B(As)={(L1).(-1),(-L1).(-L-))}- 

显然 Y(t)= 钱 。 为 了 证 明 近 似 解 Y(1) 收敛 于 方程 (5-1D) 的 精确 解 忒 (0) ， 我 们 需要 一 
些 引 理 : 

引 理 SIM Gatv, HFIS Joj E{-L0}, krok, 表示 一 1 在 多 维 指标 中 的 位 置 ， 


元 [对 可 分 解 成 积分 和 的 形式 ， 即 


P,-) 
Lpt= DN Ta, [x] HAPT, [af 


引 理 5.2 ER 5.2 下 ， 对 于 方程 (5-1) 的 解 于， 存在 非 负 常数 C, ， 使 得 


r| soplx of <c (1+ 2] of). 


引 理 5.37514 
A 2 
Gis =E ae la [x]. A > 


MF ae pr{y}, Mae{(j, jr0J,):J, E{-L01}1<Si<JJEN}. SEF TAK, 
对 于 任意 的 站 <p<s<7，c-p<A<l， 存在 常数 全 了 和 (4+)， 使 得 


a WatPa SPa Adatna- f° pv 
Go, < Aretha NPA J G° dt. 


引 理 5.42514 
k,- 
E EDA ASAN, s| 
sejo) | n=0 E h 


对 于 aEf\{v}, Mac{(i jj): j, E{-10.1},1<i<J, JEN}, G, <oo, 分 割 
(t) ={2,:0=4, St S- ST, KEN, supjey (ti -4) SA}, Hk, =max{keN:t, <T}. 对 
于 固定 的 a ， 存 在 有 限 正常 数 C, (o) ， 使 得 


t 
. tA n Gy ds A Wa = Pa =0 
Veort < 


C, (a) A f" Gds :w+p,>0 
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其 中 Ge = =| swt DA 本 
引 理 5.5 在 假设 5.1 和 假设 52 下 ，0<t<7T， 近 似 解 Y(1) 满 足下 面 的 不 等 式 


zfs of 辐 sc (a). 


证 明 : 根据 (5-32) 可 得 
k—l 
erorile BE ahenn] 


ait (tX) II" +2. en 


+> [Sale pY IY +2, lle (42%) 


ae B{A;)\ i=0 
依据 基本 不 等 式 可 得 
Esref 加 < ERP + Je, E| sup 2 A GA Ra ills 


O<s<t 


+ dF 


aca, UB(A;) 


D keri + ln (nk x| 


sup 
O<s<t 


| 


应 用 引 理 .5.3 的 结论 ， 得 到 


[splrGs (sf <n sup 


O<r<s 


2 
a A 
lal ta 


i 


2 


En? + p> D,( a) f 


Ay} 


2 


+ 25 D,(a) f E < 


sup |}, 
a€A,UA(A;) OSrss 


T 
A 


A 


[soplr (ef Is]<D A+ > Dla) a)K? (a) Í E if 


1+ sup 


O<r<s 


se 
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2 


ds (5-33) 
= 


1+ sup 


0<r<s 


+ >> Dl(a)k lo) fz F 


aCA UB(A;) 


Z| ， 依 据 (5-28) 式 ， 可 得 


2 


= E| sup |Y 


0<r<s 


4 


r +|, pÉ, A 
A a e 
the, A 
Al lA 


A 
1+£ 


2 2 


<2E| sup +2E 


0<r<s 


F, 


A 


A sup 


0<r<s 


S$ 


2 


JE FZ |+24° K? Y 


sup |Y 
A 


O<r<s 


sup 
O<r<s 


4 


s (5-34) 


路 
a 让 


<D,E|¥,| + {F 2. D(a)Ki(a)+ > 2D,(a)Kj(a)(1+A°K?) (aah 
acA\{ @eA,VB(A3) Osrs 
t 
< D,E|Y,} +D, (æ) P| 
0 O<r: 


AG sja 


< D.E\¥,{ + D, (æ)T +D, (a) Je(sup IY (A | Jas 
0 isrss 


A 


d 


2 


sup |Y 
0<r<s 


=24°K; +2(1+.A°K?)E 


将 (5-34) 代 入 (5-33)， 得 出 


2 


a8 
+ >》 D,(a)K;( cfs od 


1 
(soplr (sf Iz) <D, 区 + >》 D,(a)K;(a) at [e 
Ossst aeA,\{v} 0 Osrs. 


2 


@eA,VUf(A;) Osrss 


2 


2 


接着 ， 应 用 Gronwall 不 等 式 ， 就 可 得 到 
¥(s} |< DEl +D, )T +D,( (a) f e” (DEK) +D, (a)i 


0 


Bfsup 


0<s<t 


< (DE| +D, (a)T) e> 
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一 Ci (oj 
定理 5.1 在 假设 5.1 和 假设 5.2 F, IFAS- X (t) PAREY (t) ， 存 在 非 
负 常 数 C ， 使 得 


E| sup |X (t)—Y (t)) I<c.a. 
[sn -re |%|s 


证 明 : 根据 (5-32) 和 (5-21) 可 得 
x (0) -ro -koy F (Stell xd + ,| 
+ i, LEDI +L (XN, JX 


aeB(A) 


~% (Se AOIRA (Hs 这 | 
-ZE (See ez, ] 


a Eå, 


BS te) 


i=0 


2 


A A eX), | 


acA\fy} 


A 60 +4[.62)-1.64], 


aeB(A; 


+ (aeoo) 
+E ae [u (7) -L(Y "+t [7)- (tX )] | 


aep(A,) 


+ 5 (Enh ener + ult) -u(x ) 


CEA2 


2 


(5-35) 
对 (5-35) 式 ， 两 边 平 方 后 先 取 上 确 界 ， 再 取 条 件 期 望 得 
E su O-ro z | 
<al a+ da] (5-36) 
acA\{y} aeB(A;) aef{A,) aeB(A;) aeA, 


其 中 
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及 = 五 ae 


te[0,7] 


Soll eX) Ctl A tonX | 5} 


S* =E] su pe, [LOLLO +(X) | | 


[or] 


U, = 


A A ls} 


dal 
o -alaaf 


enti en-e) 
j ike Ta [4 (4%) -la a4 i +1, [ia (tot, J- (42% j) 


2 
fis} 


在 假设 5.1 下 ， 依 据 引 理 5.4 3 R HT, 


=E sup [$ [i (2)- E NAg +lk (tika) (到 让 3] 
< D,(a) f 中 加 1, (s,X (s))—1, (s.r (s) 5 dt 
K? (a) i a| sup |X (5) - Y(s)f 四 dt (5-37) 


在 假设 5.1 下 ， 依 据 引 理 5.4 Xt SP RET AT, 
S? < D, (a) fe{ sw 


< D, (a) K} (æ) 四 品 KO-rGj Ja (5-38) 


1, (s,X (s))-1,(s.¥(s))} Ja 


现在 根据 引 理 5.4， 估 计 Ue， 


or sma fafagh xo ea 


<D,(a) 8°? (a) f 


1+ E| suplX (s |%lla 
{spot 


<C,(a)A’ (5-39) 
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KPG, (a)= D, (a) Ki (a)T(1+C, (+z of). 
在 假设 5.1 下 ， 依 据 引 理 5.4 Xt Or 进行 估计 ， 


0° <C, (aa mp ls, ¥(s))-1 e) 中 
< 2C, (a) K} (a) as AN dt +2C, (a) K? (a) sla Y(s) “Hy 对 
< <TC, (a) K} (a) A’ +2C, (a) K} conf m Y(s)- Tey r \ (5-40) 


2 
mes ¥(8)-Fa s) 4re[t,.t.,) > K5 5.5 的 结论 代入 
sel0, 7 
(5-34) R40 E sup jY FAR. BAA GS-28)K. (5-31)K. BW 5.1、 假 设 5.2、 引 理 


5.4 和 引 理 5.5 可 知 


æo m ka hg s 


D lett E Jello, wl D elel SJN 


acA \{v} aep(A;) 


vel 


<D, 
aeA\\(y} (x 


“oat 


< s sup 
seft, .t] 


a (tas, 


Y, Ai | |s }_ > E © E| spz sup |Z 
aeA,UP(A;) sefi J 


af 


"| %)as+D, Jeh (x) 


efir] 


< Bod felito 
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<D,A*(1+2{ i F))+Ba(1+# (ye 列 
<B,(a)A aiii 
将 (5-41) 带 入 (5-40) 得 


O? <217¢, (a)K3(a)A*-+2C, (a) Ba (a) K? (a)? 
<6,(0)A? (5-42) 


在 假设 5.1 和 假设 5.2 下 ， 利 用 式 (5-28)、 引 理 5.4、 引 理 5.5 对 Ze 进行 估计 ， 


2 
sup 5 | 
| 罗 T] 


Sf (ts) (下放 (天 让 


i=0 


<D, “ai (zo). BSA |a 
< D, (a) K$ (a) i Éa) ig sa 
< D, (a) K? oj Aa) sha 


2 


Y 


< D, (a) K? (a) K2A? fi 站 中 


<C,(a)A’ (5-43) 
FEC, (a) = D, (a) Ki (a) K? (1+ D,(a))T 。 
将 (5-37)-(5-39) 和 (5-42)-(5-43) 代 入 (5-36) 可 得 


E apir (t)- v(t 加 <C,(a) 


(D, (a) +D,(a Ka) [sap to)-r (sf |G] ar+6 (a) A* +6, (a)A +6, (a)A 
sE an (s)—Y( (sf A 


运用 Gronwall 不 等 式 可 得 


华东 理工 大 学 硕士 学 位 论文 第 61 页 


E| sup |X (t)—¥ (t |%|<Cd?, 


rel0,7] 
其 中 C=C, (a)r. 


定理 5.1 表明 在 均 方 意义 下 离散 解 了 (1) 以 速度 1 KAFEI X (1). XRS 
出 了 欧 拉 算 法 的 一 阶 显 式 ， 即 
Vin =Ye+f (bh Fe)Ate(h, Ve )AW, +h(h, Ye )AN, 


+58 (tot). ((am,) -AN,) 


+ 人 (下 +h(to Yi ))-a(t Y) (AN7 -AN, } (5-44) 
+(g (hs 下 +h(t,, Y))-2(4.¥"))AMAN, 


+(e (tY. (A (4%) +1)-2(t% salt) È (Mew -™ ) 


Na +1 ‘t 
下 面 用 一 个 例子 来 比较 本 文 的 Split-step 算法 (5-28)、(5-29) 与 欧 拉 算法 (5-44)。 令 
T=1,X,==10, 和 =50。 方 程 (5-1) 的 系数 分 别 为 了 (1,x)= 2x,g(t,x)= Sx h(i) = -i ， 
方程 (5-1) 的 解析 解 为 


oy 15 1 
X(t)}=10|— —t+—Wl. 
人 a apl my | 


5.1 两 种 算法 得 到 的 近似 解 和 精确 解 的 路 径 
Fig. 5.1 Exact Solution, Solution of the order 1 Euler-Maruyama scheme and of the order 1 Split-step 
scheme 
文中 所 给 出 的 近似 解 和 误差 是 100 次 计算 结果 的 均值 。 图 5.1 给 出 A=0.01 的 
Split-step 算法 的 一 阶 近 似 解 和 Euler-Maruyama 算法 一 阶 近 似 解 及 精确 解 , 从 图 5.1 可 以 
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看 出 Split-step 算法 的 数值 更 近似 于 精确 解 . 图 5.2 显示 的 是 Split-step 算法 的 一 阶 近似 解 
和 Euler-Maruyama 算法 一 阶 近似 解 与 精确 解 之 间 的 误差 , 由 图 5.2 可 知 ，Split-step 算法 
的 一 阶 近似 解 与 精确 解 之 间 的 误差 更 小 ,此 例 说 明 Split-step 算法 要 优 于 Euler-Maruyama 
算法 。 


wa F 
Ah 


4 


5.2 离散 解 与 精确 解 取 对 数 后 的 误差 
Fig. 5.2 Point mean square errors of the Euler scheme and of the Split-step scheme in a logarithmic scale 
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61 总 结 


本 文 给 出 了 三 种 求解 带 跳 的 随机 微分 方程 数值 解 的 算法 。 针 对 一 类 带 跳 的 随机 微分 
方程 ， 本 文 给 出 了 基于 Euler-Maruyama 法 的 Split-Step 算法 ， 证 明了 当 方程 的 系数 满足 
Lipschitz 条 件 和 线性 增长 条 件 , 且 时 滞 函 数 满足 一 定 的 条 件 时 ，Split-Step 算法 的 收敛 速 


Ran yaz ， 并 通过 两 个 数值 实例 验证 了 算法 的 有 效 性 。 随 后 ， 对 不 含有 时 灌 项 的 
. Split-Step 算法 做 加 权 得 到 了 SoM, 证 明了 当 方 程 的 系数 满足 Lipschitz 条 件 和 线性 增 
长 条 件 时 ， S59 算法 近似 解 的 收 人 速率 是 了 。 然后 , 将 SS9 算 法 推广 到 了 带 时 滞 项 的 情况 


Fo 证 明了 带 跳 的 随机 延迟 微分 方程 近似 解 的 收敛 速率 是 。 在 最 后 一 部 分 ， 本 文 应 用 


It6—Taylor 展开 式 对 Split-Step 算法 做 二 阶 展开 ， 得 到 16 一 Taylor 展开 的 近似 解 ， 并 证 
明了 所 得 的 近似 解 的 收敛 速率 是 1。 
6.2 展望 


在 文章 的 最 后 , 我 们 提出 关于 随机 微分 方程 数值 解 方 面 的 工作 今后 可 以 做 如 下 展开 : 

(1) 我 们 的 方程 所 考虑 的 泊 松 跳 是 一 个 泊 松 过 程 , 但 实际 上 更 一 般 的 可 以 考虑 将 跳 过 
程 推 广 到 泊 松 跳 测度 。 

人 2) 我 们 在 证 明 Split-Step 算法 一 阶 近似 解 的 收敛 速 率 时 ， 并 没有 考虑 时 滞 项 ， 因 此 
随机 微分 方程 中 引入 时 滞 项 也 是 可 以 继续 探讨 的 工作 之 一 。 
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